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Preface 


This  independent  study  began  as  an  effort  to  perform  a 
follow-up  analysis  of  another  thesis  on  airborne  penetration 
of  radioactive  clouds.  The  original  thesis  was  done  by  B.E. 
Hickman  (Ref  10)  and  performed  a  worst  case  analysis  by  compu¬ 
ting  the  doses  aircrew  members  would  receive  when  flying 
through  a  descending  fallout  cloud  generated  by  a  nuclear 
surface  burst.  Military  planners  are  still  interested  in  this 
problem,  due  to  the  proposed  dense  packing  of  the  MX  missiles. 
An  attack  on  these  silos,  housing  the  MX  missiles,  by  megaton 
or  greater  sized  weapons  would  present  a  serious  problem  for 
strategic  aircraft  carrying  out  their  wartime  missions.  Since 
these  aircraft  would  be  airborne  for  extended  periods  of  time, 
penetration  of  radioactive  clouds  would  be  probable. 

Presented  within  are  my  efforts  on  extending  and  improving 
the  previously  developed  computer  code.  This  code  calculates 
the  ionizing  dose  rate  caused  by  the  portion  of  the  radioactive 
cloud  ingested  by  an  aircraft  as  a  function  of  time.  It  also 
computes  the  dose  rate  caused  by  radiation  external  to  the 
aircraft.  The  external  radiation  will  be  hereafter  referred 
to  as  sky-shine  radiation.  Radiation  doses  from  a  multiple 
burst  scenario  are  also  computed  by  an  extension  of  this  code. 

I  am  especially  thankful  to  Dr.  Charles  J.  Bridgman  for 
his  patience  and  directi cn  during  this  research.  I  am  also 
very  grateful  for  the  patience  and  love  my  wife,  Karen,  showed 
during  the  period  of  this  research  work. 
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This  report  evaluates  the  threat  to  aircrew  members 
when  their  aircraft  approaches  and  subsequently  penetrates  a 
descending  radioactive  cloud  generated  by  a  nuclear  weapon 
surface  burst.  — _ 

The  computer  code  developed,,  during  this  study,  is  a 
revision  of  the  code  developed  by  Hickman.  \  The  re-development 
of  Hickman's  program  consists  of  a  remodeling  of  the  computa¬ 
tional  methods  for  sky-shine  dose  and  cloud  model.  The  code 
also  computes  the  ionizing  dose  rate  an  air  crew  member 
receives  when  flying  through  the  radioactive  cloud  as  a  func¬ 
tion  of  time. ^The  cloud  model  developed  is  patterned  after 
the  AFIT  (Fallout  Smearing)  code.  A  comparison  is  made  between 


the  activities'7 and  dos**s  received  between  500  and  12,000  meters 
altitude  at  increments  of  500  meters. 

^The  code  computes  the  doses  by  considering  the  cloud 
size,  the  aircraft's  transit  time,  the  ingestion  rate  of 
radioactive  particles,  the  aircraft's  distance  to  the  burst, 
and  the  aircraft's  altitude.  A  simple  extension  of  the  com¬ 
puter  code  computes  the  dose  received  from  multiple  bursts. 

,Lirtr,  <."t 

The  results  show  that  at  9500  meters  ( approximately  31,000 

A  aJjtrvX. 

feet) ,  the  total  dose  to  each  aircrew  member  is  approximately 

A 

5  rem  after  flying  through  the  cloud  1  hour  after  cloud  sta¬ 
bilization.  The  multiple  burst  dose  is  approximately  204  rem 
under  the  same  conditions  as  the  single  burst  case.  Both  the 
single  and  multiple  burst  case  use  a  mission  completion  time 
of  8  hours  after  entering  the  cloud. 


AIRBORNE  PENETRATION  OF 


RADIOACTIVE  CLOUDS 

I.  Introduction 


Background 

The  proposed  basing  of  the  NX  missile  system,  whether  by 
dense  packing  or  a  similar  scheme,  would  definitely  be  a  prime 
target  for  an  enemy  first-strike  attack.  The  numerous  radio¬ 
active  dust  clouds  generated  by  such  a  strike,  aided  by  the 
prevailing  winds,  would  cause  large  areas  of  the  continental 
United  States  to  be  exposed  to  radioactive  fallout.  Airborne 
aircraft  would  also  be  exposed  to  these  clouds  when  either 
entering  or  approaching  them.  Strategic  aircraft,  performing 
their  wartime  mission,  would  be  especially  susceptible  to 
exposure.  The  Strategic  Air  Command  has  expressed  concern 
about  the  amount  of  fallout  that  Airborne  Command  Post  per¬ 
sonnel  would  be  exposed  to  when  flying  through  these  radio¬ 
active  clouds  (Ref  10) .  The  specific  dose  rates  generated  by 
this  study  provide  information  to  help  redirect  an  aircraft’s 
flight  path  if  location,  time  of  burst,  approximate  yield  and 
wind  profiles  are  known.  By  knowing  these  parameters,  the 
aircraft's  altitude  could  be  decreased  in  order  to  reduce  the 
dose  caused  by  gamma-ray  radiation  to  the  aircrew. 

Aircraft  penetration  of  radioactive  clouds  is  hazardous 
in  at  least  four  ways.  First,  the  aircrew  is  exposed  to 
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ionising  radiation  through  the  aircraft's  ventilation  system 
and  through  the  aircraft's  skin.  Secondly,  the  aircrew  may 
injest  the  radioactive  particles.  Thirdly,  electronic  equip¬ 
ment  could  malfunction  if  the  ionising  dose  rate  is  high 
enough  (Ref  9) .  Fourth,  if  the  dust  particles  are  large  and 
numerous  enough  the  aircraft's  engines  could  be  physically 
damaged  by  mass  ingestion  of  the  dust  particles.  This  study 
focuses  only  on  the  first  hazard.  Note  that  the  second 
hasard  could  be  all  but  eliminated  by  having  the  aircrew  put 
on  their  oxygen  masks. 

Problem 

No  data  on  previous  flights  through  radioactive  clouds 
could  be  found.  Therefore,  comparisons  are  not  made  against 
actual  data.  Comparisons  are  made  against  the  data  generated 
by  Hickman's  (Ref  10)  computer  code  and  were  found  to  be 
approximately  20  percent  lower  than  his  for  a  single  meagaton 
else  burst.  The  parameters  used  in  this  study  were  for  an 
SC-135  aircraft.  This  is  the  aircraft  used  as  an  airborne 
command  post.  These  parameters  are  the  same  for  a  KC-135 
(refueling)  aircraft.  Parameters  for  a  B-52G  could  also  be 
used  in  the  computer  program  (Ref  17) .  Although,  because  of 
time  limitations,  the  author  did  not  compute  the  dose  to  a 
B-52G  aircrew  member. 

The  computer  program,  developed  during  this  study,  com¬ 
putes  particle  size  distribution  as  a  function  of  yield, 
altitude  and  fall  time  for  a  specific  time.  With  the  particle 
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fiiae  distribution  known,  from  sea  .level  on  up,  the  amount  of 
activity  at  any  altitude  can  be  computed.  Cabin  dose,  caused 
by  the  ingestion  of  particles  at  the  aircraft's  altitude,  and 
sky-shine  dose  can  then  be  computed.  Total  aircrew  dose,  due 
to  gamma  radiation,  will  then  be  known. 

Scope 

This  study  models  the  radioactive  cloud  similarly  to 
Hickman's  approach.  This  computer  model  computes  both  the 
activity  that  enters  the  aircraft  and  the  radiation  directed 
through  the  external  skin  at  specific  times  and  altitudes. 
There  are  four  major  changes  compared  to  Hickman's  model. 
First,  the  sky-shine  dose  computational  method  is  changed  as 
to  not  incorporate  Hickman's  approximation,  but  rather  a  dif¬ 
ferent  approximation.  Hickman  computed  sky-shine  dose  by 
assuming  that  the  aircraft  was  positioned  at  the  center  of 
the  radioactive  cloud  for  the  time  it  took  the  aircraft  to 
fly  through  the  cloud.  This  study  computes  the  sky- shine 
dose  by  integrating  the  normalized  horizontal  distribution 
function,  f(x,t),  from  minus  to  plus  infinity.  Second,  par¬ 
ticle  size  distribution  centers  are  net  all  set  to  the  cloud' 
center  altitude  at  stabilization  time.  The  distribution 
centers  are  instead  set  to  individual  particle  altitudes  at 
stabilization  time.  Third,  a  shielding  approximation  was 
used  in  computing  the  sky-shine  dose.  Fourth,  integrati.w 
from  minus  to  plus  infinity  of  the  f (x, t)  function  was  also 
used  to  compute  the  activity  caused  by  the  radioactive 
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particles  that  enter  the  cabin.  Hickman,  on  the  other  hand, 
used  two  times  the  standard  deviation  of  the  f(x,t)  function 
as  the  lateral  limits  of  the  radioactive  cloud.  His  program 
then  computed  the  amount  of  radioactivity  that  entered  the 
aircraft  during  the  time  it  took  the  aircraft  to  traverse 
the  radioactive  cloud.  Also,  for  internal  activity,  the 
numerical  integration,  used  to  compute  the  dose  rate  at  the 
center  of  the  cabin,  used  a  different  value  for  cabin  radius 
and  consequently  a  different  value  for  this  integral  was 
calculated  and  used  in  the  computer  code.  These  changes  are 
discussed  in  greater  detail  later  in  the  text.  The  results 
for  cabin  activity  still  constitute  a  worse-case  approach. 

This  is  because  all  of  the  particles  that  are  ingested  into 
the  cabin  are  assumed  to  remain  in  the  cabin  for  the  duration 
of  the  mission.  The  accumulation  of  radioactive  particles  in 
the  cabin  is  due  to  cabin  pressurization,  air  conditioning 
and  electronics  cooling.  Only  the  mass  flow  rate  from  cabin 
pressurization  and  air  conditioning  is  used  in  this  study. 

A  different  compressor  with  a  different  mass  flow  rate  is 
used  for  electronics  cooling  (Ref  19) .  The  aircrew  are  also 
exposed  to  radiation  by  direct  contact  and  inhalation,  although 
only  tissue  dose  from  the  external  gamma  radiation  is  addressed 
in  this  report. 

Assumptions 

There  are  several  explicit  assumptions  made  in  this 
report.  They  are: 
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1.  The  rky -shine  dose  is  integrated  between  the  times 
of  plus  infinity  a ad  minus  infinity 

2.  The  initial  conditions  for  the  stabilized  cloud  as 
treated  by  Bridgman  and  Bigelow  (Ref  1)  and  remodeled  by 
Hopkins  (Ref  2)  are  applicable 

3.  The  particle  size  distribution  at  cloud  stabiliza¬ 
tion  time  as  treated  by  Bridgman  and  Hickman  (Ref  2)  is 
applicable 

4.  There  is  no  self  shielding  of  gamma  radiation  inside 
the  cabin 

5.  There  is  no  adherence  of  radioactive  particles  to 
the  exterior  of  the  aircraft 

6.  Megaton  size  thermo- nuclear  weapons  will  be  used 
against  the  hardened  missile  silos. 

These  assumptions  are  also  discussed  in  more  detail  later  in 
the  text. 

Approach 

The  mathematical  development  of  the  computer  model  and 
a  summary  of  the  results  for  a  single  and  multiple  burst 
scenario  in  terms  of  activity  density  in  Curies  per  cubic 
meter  versus  altitude  at  various  times  are  presented  in 
Chapter  XI.  The  mathematical  development  for  the  external 
dose,  caused  by  both  cabin  activity  and  sky-shine,  is  pre¬ 
sented  in  Chapter  III.  The  results  for  a  single,  1-megaton 
ground  burst  are  then  presented  in  tabular  form.  These 

5 


V> 


*«***’• 

-  ■  ■  . 


tables  include  the  doses  received  and  the  particle  contri¬ 
buting  the  most  activity  at  the  specified  altitude.  The 
mathematical  development  for  incorporation  of  multiple  burst 
doses  into  the  computer  code  is  presented  in  Chapter  IV. 

The  results  are  then  presented  in  tabular  form  for  a  specific 
multiple  burst  scenario.  Conclusions  and  recommendations 
are  presented  in  Chapter  V. 


■\  *-  .v.v . 


models  are  WSEG-10  and  DELFIC  (Ref  20) .  WSEG-10  uses  empiri¬ 
cal  functions  rather  than  numerical  analyses  and  computes 
fallout  dose  rates  in  seconds  or  less  on  modern  computers. 


WSEG-10  is  normally  used  in  operational  type  studies.  DELFIC 
is  used  for  research  work  and  as  a  comparison  standard  because 
of  its  excellent  predictive  capabilities.  Unfortunately, 


DELFIC  is  cumbersome  to  use  and  slow  to  compute.  On  the 
other  hand,  WSEG-10  fails  to  account  for  fractionation,  to 
allow  for  activity  particle  size  variations  in  the  fallout 
distribution  and  to  use  realistic  settling  rates. 

Patrick  (Ref  16)  and  Patrick  et  al  (Ref  17)  investigated 
aircraft  penetration  of  radioactive  clouds  using  a  simplified 
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model.  This  model  was  perhaps  too  simple,  since  it  treated 
radiation  variations,  fractionation  and  particle  size  dis¬ 
tribution  rather  arbitrarily. 


This  study,  like  Hickman's,  uses  the  fallout  prediction 
model  developed  by  Bridgman  and  Bigelow  (Ref  10) .  This  is 
the  AFIT  model  and  produces  results  which  compare  favorably 
with  DELFIC  (Ref  1) .  Since  this  model  accounts  for  fraction¬ 
ation,  activity  size  variations,  uses  realistic  settling 
rates  and  computes  in  seconds  or  less  it  is  used  in  the 
development  of  the  model  in  this  study.  Since  only  the 
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airborne  cloud  is  used  in  this  study,  the  g(t)  function, 
which  is  the  normalized  rate  of  arrival  of  activity  on  the 
ground,  was  not  applicable. 

Theory 

The  DELFIC  model  is  considered  to  be  the  standard  pre¬ 
dictive  model.  The  cloud  model  presented  by  Bridgman  and 
Bigelow,  as  modified  by  Hopkins  (Ref  2) ,  is  used  in  this 
study.  Activity  size  distributions  are  used  which  do  account 
for  fractionation.  These  distributions  are  divided  into  100 
discrete  size  groups.  The  size  groups  used  in  this  study 
are  from  the  DELFIC  default  parameter  (Ref  1) . 

The  cloud  produced  by  a  nuclear  burst  consists  of  hot 
gases  and  vapors  which  rise  immediately  from  the  detonation 
point.  The  cloud  is  cooled  by  radiation  and  convection  until 
its  contents  reach  the  local  ambient  atmospheric  temperature 
at  which  point  the  rise  ceases.  The  cloud  stabilizes  in  4  to 
10  minutes  at  a  height  and  radius  determined  by  the  weapon ' s 
yield  (Ref  9) .  Suspended  solid  particles  are  formed  during 
the  cooling  phase  and  before  cloud  stabilization.  These  par¬ 
ticles  are  approximately  spherical  in  shape  with  sizes  ranging 
from  less  than  a  micron  to  several  millimeters  for  surface 
bursts. 

The  radioactive  cloud  is  modeled  as  an  upright  volume 
which  approaches  a  right  circular  cylinder.  This  assumes  no 
variation  in  wind  translation  from  the  top  to  the  bottom  of 
the  cloud.  The  radioactivity  contained  within  the  cloud  is 


assumed  to  be  normally  distributed  in  the  horizontal  direc¬ 
tion  and  individual  size  groups  are  assumed  normally  distri¬ 
buted  in  the  vertical  direction.  Although,  these  assumptions 
are  almost  never  true,  they  don't  really  matter  for  two  rea¬ 
sons.  First,  the  sky-shine  dose  is  determined  by  the  particle 
size  activity  groups  that  are  near  the  aircraft.  Distant 
groups  are  unimportant  due  to  the  exponential  attenuation  of 
air.  Second,  for  the  cabin  activity,  the  aircraft  flies 
through  the  whole  cloud  at  a  constant  altitude  sweeping  out 
all  of  the  particles  that  are  at  that  altitude  at  that  time. 
Figure  1  illustrates  this  model. 

The  radioactivity  per  unit  volume  of  outside  air  is  a 
function  of  particle  radius,  r,  and  time,  t.  The  activity 
density  is 

00 

A"'(x,y,z,t)  =/  A*  •  '  (x,y,?.,r,  t)  dr  [Ci/m3]  (1) 

o 

Variation  of  the  specific  radiation  density  term,  Ar' 1 ',  in 
any  of  the  three  spatial  directions  is  assumed  to  be  inde¬ 
pendent  of  the  other  two.  It  is  also  assumed  that  variation 
in  the  two  horizontal  directions  is  independent  of  particle 
size.  Therefore,  Eq  (1)  becomes 

Ar '  '  '  (x,y,  z,  t)  =  f(x,  t)f(y,t)  /  A  ' (z, r , t) dr [Ci/m3]  (2) 

o  l 

where 
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gure  1.  Fallout  Cloud  Model 
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Where  ox  and  cy  are  defined  later  in  this  section.  By  re¬ 
placing  the  integral  in  Eq  (2)  with  a  summation  over  the  100 
discrete  particle  size  groups,  where  each  group  is  assumed  to 
consist  of  a  mono-sized  particle,  r^,  which  contains  one  per¬ 
cent  of  the  total  activity  one  obtains 


/°°  Ar '  (z,  r,  t)  dr  [Ci/mJ  =  £  A<  (r,  t)  f i  (z,  t) 
o  L  i=l  A 


Where 


fA(z,t)  = 


/2ir 


fi?)’  [;] 


(5) 


(6) 


In  equations  (3)  through  (6)  x  is  the  downwind  distance, 
y  is  the  crosswind  distance,  is  the  vertical  height  for 
each  particle  sxze  group,  z  is  the  aircraft  altitude,  here¬ 
after  za,  and  the  sigmas  are  the  time  varying  standard  devia¬ 
tions.  Time  of  arrival,  ta,  was  used  in  the  program  instead 
of  true  time,  t.  This  is  an  approximation  since  the  time 
does  vary  as  the  aircraft  traverses  the  cloud.  But  since 
this  time  is  from  one  to  two  minutes,  ta  approximately  equals 
t  during  aircraft  passage  of  the  cloud.  In  WSEG-10,  v  is  a 
constant  wind  which  rotates  uniformly  with  altitude.  In  this 
report  the  aircraft  velocity,  VAC,  is  used  instead  of  vx. 
Therefore,  since  the  aircraft  and  the  cloud  are  in  the  same 


11 


medium,  VAC  is  the  true  airspeed,  not  groundspeed,  of  the 

aircraft. 

Eq  (4)  can  now  be  constructed  as 

A'"(x,y,z,t)  =  f  (x,  t)f  (y,  t)  12°  A.  (r ,  t)  f1  (z,  t)  (7) 

i-1  1 

The  activity  A^(r,t)  is  equal  to  one  percent  of  the  total 
cloud  activity.  This  was  accomplished  by  evaluating  the 
cumulative  activity  size  distributions  after  running  the 
computer  code  DELFIC  with  its  default  particle  size  distri¬ 
bution  (Ref  2) .  Other  activity  size  distributions  could  be 
constructed  by  the  method  described  by  Bridgman  and  Bigelow 
(Ref  1) .  The  mean  radii  of  the  100  particle  size  groups  are 
given  in  Table  I.  The  time  variation  of  the  activity  follows 
the  Way-Wigner  approximation,  A(t)  =  Aj^  x  f1*2,  where  A(t) 
is  the  cloud's  activity  at  time  t  and  A^  is  the  activity  at 
unit  time.  530  gamma  megacuries  per  kiloton  of  fission  yield 
was  used  for  the  activity  at  unit  time  (Ref  9) . 

Computation  of  both  the  sky-shine  and  cabin  ingestion 
doses  performed  by  integrating  the  normalized  horizontal 
distribution  function,  f(x,t),  from  minus  :o  plus  infinity. 

This  insures  that  all  of  the  particle  size  groups,  located 
at  the  aircraft's  altitude,  contribute  to  the  dose  received. 
Therefore,  this  study  does  not  use  a  worst  case  approach,  as 
Hickman's  did.  The  methods  used  to  confute  sky-shine  and 
cabin  ingestion  doses  are  detailed  in  their  respective  sections. 
The  assumptions,  which  allow  the  above  methods  to  be  used. 
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TABLE  I 

Mean  Radii,  in  Microns,  of  the 
100  Equal  Activity  Groups 


m 


.48 

.92 

1.29 

1.64 

1.99 

2.35 

2.70 

3.07 

3.45 

3.83 

4.23 

4.64 

5.06 

5.49 

5.94 

6.39 

6.87 

7.36 

7.87 

8.39 

8.93 

9.50 

10.1 

10.7 

11.3 

11.9 

12.6 

13.3 

14.0 

14.7 

15.4 

16.2 

17.0 

17.9 

18.7 

19.5 

20.5 

21.5 

22.5 

23.5 

24.6 

25.7 

26.8 

28.0 

29.3 

30.5 

31.9 

33.2 

34.7 

36.2 

37.7 

39.3 

41.0 

42.8 

44.6 

46.5 

48.4 

50.5 

52.7 

55.0 

57.3 

59.8 

62.5 

65.3 

68.2 

71.3 

74.5 

77.7 

81.4 

85.2 

89.2 

93.5 

98.0 

103. 

108. 

113. 

119. 

126. 

132. 

140. 

148. 

156. 

166. 

176. 

188. 

201. 

215. 

231. 

249. 

270. 

294. 

323. 

357. 

398. 

450. 

519. 

615. 

762. 

10  3  3. 

1873. 

r* 

,\ 
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are  also  defined  in  these  sections.  Methods  for  computing 
sigma  x  and  y  have  been  developed  for  both  the  WSEG-10  (Ref 
18)  and  APIT  (Ref  1)  models.  Although,  since  the  g(t)  func¬ 
tion  was  uot  used  sigma  x  was  never  calculated  in  these  two 
reports.  It  is  assumed  that  sigma  x  will  behave  as  sigma  y, 
except  without  the  torroidal  growth  term.  This  term  is  not 
needed  for  sigma  x  since  a  constant  wind  in  the  x  direction 
is  assumed.  The  equations  for  sigma  x  and  y  are 

°v(ta>  =  °o  I1  +  +  [o2*s,t a1  (8> 

TC 


where 


°^(ta)  *  a*  [1  +  8t^] 

TC 


(9) 


eg  *  exp [0 . 7+Jln  (YLD)  /3  -  3. 25/ (4+ Un  (YLD) +5. 4] )  J  (10) 


and 

-  (HC)  2 

Tc  =  [12fiC  -  (2.5H£)  ]•  [1-0.5  e  625  1  (H) 

60  60 


where 


HC  *  44.0  +  6 . 1  •  Jin  (YLD)  -  0 . 205  •  [  Jin  (YLD)  +  2.42] 

•  |  Jin  (YLD)  +2 . 42  |  (12) 

The  time  constant  used  in  WSEG-10,  Tc,  and  ta  are  in 
hours  and  the  standard  deviations  are  in  meters.  HC  is  the 
center  height  of  the  cloud  and  is  also  from  WSEG-10.  HC  is 
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in  kilofeet  and  YLO  is  in  megatons. 


S  is  the  wind  shear  and 


is 


(13) 


Shear  is  assumed  to  have  a  constant  value  of  1  in  this  study. 
Sigma  z  will  be  defined  later  in  the  text.  The  time,  ta, 
will  always  be  less  than  or  equal  to  3  hours,  since  this  is 
when  the  torroidal  growth  term  is  assumed  to  die  out  (Ref  18) . 

F (z, t)  is  evaluated  by  initially  assuming  a  vertical 
distribution  of  the  cloud  at  stabilization  time  for  each 
activity  size  group  and  then  following  the  distribution  as 
it  fell  to  the  ground  (Ref  2) .  The  stabilized  time  vertical 
distributions  are  assumed  to  be  normal  functions  with  the 
lower  distribution  centers  at  the  lower  altitudes  for  the 
larger  particles. 

The  stabilized  altitude  for  the  distribution  centers 
for  each  activity  size  group,  z^,.is  calculated  from  the 
following  formula  from  Hopkins  (Ref  2) 


c2°i 


(14) 


where 

tn^  «  7.889  +  0.34-lnY  +  0 . 001226*  UnY)  2 

-  0.0U5227-  UnY)  3  +  0.000417- UnY)  4  (15) 

and 

£nC2  =  1.574  -  0.01197-inY  +  0.03636- UnY) 2 

-  0.0041- UnY)  3  +  0.0001965-  UnY)  * 


(16) 


The  initial  positions  of  some  selected  groups  are  show  con¬ 
ceptually  in  Figure  2.  is  the  particle's  diameter  in 
microns  and  Y  is  the  yield  in  kilotons. 

Hopkins  developed  Eq  (14)  by  executing  the  DELFIC  cloud 
rise  model  using  the  default  particle  size  parameter  for 
several  yields  and  then  calculating  the  altitude  of  each 
particle  si2e  group  at  stabilization  time.  He  then  empiri¬ 
cally  fit  the  vertical  mean  for  each  size.  Figure  3  illus¬ 
trates  this  method. 

The  vertical  distribution  for  each  group  at  stabiliza¬ 
tion  time  was  assumed  to  be  distributed  normally  with  a 
standard  deviation  of 

az  *  0.13  Z*  (17) 

This  is  the  method,  used  by  WKEG-10,  to  compute  vertical 
standard  deviation;  except  all  particles  start  at  the  same 
initial  altitude  in  WSEG-10.  Even  though  DELFIC  does  not 
predict  normal  distributions  for  each  group  the  results  from 
this  assumption  are  very  close  to  DELFIC 's  (Ref  2). 

From  the  initial  position  of  each  activity  size  group, 
each  of  the  100  particle  size  distributions  were  allowed  to 
fail  toward  the  ground  until  the  aircraft  intercepted  the 
cloud.  Fall  velocities  were  computed  by  the  method  of 
McDonald  (Ref  15),  using  Davies'  (Ref  5)  polynomials  with  a 
U.S.  standard  atmosphere  divided  into  a  group  of  layers  of 
constant  density  and  viscosity.  The  distance  a  particle 
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falls  In  six  minutes  was  selected  as  the  layer  thickness. 

This  time  was  selected  by  running  the  fall  dynamics  section 
of  this  study's  computer  code  and  computing  fall  velocities 
for  particles  from  1/2  to  1800  microns  in  radius  at  altitudes 
ranging  from  1000  to  10,000  meters.  Particles  of  less  than 
100  microns  in  radius  fell,  in  six  minutes,  from  slightly 
over  100  meters  to  much  less  than  100  meters,  depending  on 
the  altitude.  For  example,  a  62.5  micron  radius  particle  at 
8000  meters  altitude  fell  at  the  rate  of  1.02  meters  per 
second  and  a  1000  micron  radius  particle,  at  the  same  alti¬ 
tude,  fell  at  the  rate  of  16.9  meters  per  second. 

McDonald  showed  that  the  Reynold's  number  squared  times 
the  drag  coefficient,  Cd,  of  spheres  falling  through  a  viscous 
medium  is  a  function  of  particle  radius,  r±,  air  density, 
pair*  and  dynamic  viscosity,  n*  This  equates  to 

R2Cd  =  32-pair-ppart-g-|i  (18) 

k 

where  Pp^rt  *s  the  particle  density  taken  as  2600  and  g 
is  the  gravitational  constant.  Davies  showed  that  the 
Reynolds  number  could  be  found  by  the  following  equations 

R  =  P2Cd/24  -  2.3363  x  10"4« (R2Cd) 2  +  2.0154  x  10“6 
•  (R2Cd)3  _  6.9105  x  10~9« (R2Cd) 4 

for  R2Cd<140 


log (R)  =  -1.29536  +  0 . 986 *  log (R2Cfl)  -  0.046677 


[log(R2Cd)]2  +  0.0011235* [log(R2Cd) ] 


for  3<R< 10 , 000 


The  local  velocity  is  found  from 


V  - 

1  2pair  ri 


Finally,  by  taking  the  velocity  of  the  particle  times  the 
fall  time  of  each  increment  the  distribution  center  of  each 
group  can  be  found  by 


Z.  (t)  =  Z^  -  E  Vi  At- 

i  -  s  j=i  i 


where  Z.^(t)  is  the  distribution  center  at  time  t  of  particle 
size  group  i,  Z  is  the  distribution  center  at  cloud  stabili- 
zation,  v|  is  the  velocity  of  particle  size  group  i  during 
one  increment  and  t  is  the  fall  time  for  that  increment. 

The  activity  density  for  a  single  burst  measured  at  the 
aircraft's  altitude,  Z_,  is 

a 


A' ' ' (x,y, z_, t)  =  f(x,t)f(y,t)  £ 


iQ0  0 . 1A ( t ] 


=1  /2tt  a; 


_k/z-za\2 


V— 
V  °i 


where  A(t)  =  — ^ 

100 

and  A*  =  530  x  106*3./  x  1010*YLD*FF 
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i  YLD  is  the  weapon's  yield  in  kilotons  and  FF  is  the  fission 

i*i*  * 

fraction.  This  activity  is  the  unit  time  reference  activity 

i 

in  disintegrations  per  second.  Therefore,  the  Way-Wigner 
|  decay  law  must  be  used  to  obtain  real  time  activity.  Eq  (22) 

|  is  the  activity  per  cubic  meter. 

I 

Cabin  dose,  due  to  ingestion  of  radioactive  particles, 
and  sky-shine  dose  will  be  covered  in  Chapter  III.  Multiple 
bursts  will  be  covered  in  Chapter  IV. 


Ill .  Cabin  Analysis 


The  parameters  used  in  this  study  are  for  an  EC-135,  an 
Airborne  Command  Post  aircraft.  These  parameters  would  be 
the  same  for  a  KC-135,  a  refueling  aircraft.  The  EC-135 
cruises  at  231.5  m/sec  and  has  a  skin  thickness  of  0.063 
inches  (Ref  17) .  The  EC-135  also  has  an  internal  volume  of 
246.357  cubic  meters  (Ref  14)  and  a  mass  flow  rate  of  outside 
air  into  its  air  conditioning/pressurization  system  of  approx¬ 
imately  150  lbm/min  (Ref  13).  The  aircrew  will  be  exposed  to 
radiation,  from  both  cabin  ingested  particles  and  from  ex¬ 
ternal  radiation  at  the  aircraft's  altitude. 

Ground-Shine 

Ground-shine,  cloud  radioactivity  that  has  fallen  to 
the  ground,  was  shown  by  Hickman  (Ref  10)  to  be  insignificant. 
This  is  true  only  if  the  aircraft  is  flying  a  few  gamma-ray 
mean  free  paths  (mfp)  above  the  ground.  A  one  MeV  gamma- ray 
has  a  mfp  of  approximately  120  meters  at  sea  level  (Ref  6) . 
Hickman  found  the  dose  rate  to  an  aircraft  flying  at  305 
meters  above  the  ground  and  through  the  cloud  to  be  equal  to 
the  ground  activity  times  1  x  10“^.  The  author  re-calculated 
Hickman's  equations  and  found  the  values  to  be  correct. 
Therefore,  ground-shine  will  be  omitted  in  the  final  analysis. 

Sky-Shine 


Approximately  6  minutes  after  a  nuclear  ground  burst 


'» 1  ^7  ^TT^.'T  fLT"*^>l  ^  r.  *"r*  ^  * 


the  radioactive  cloud  stabilizes.  The  emission  of  gamma 
radiation,  from  the  fission  products  in  this  cloud,  is  the 
only  radiation  source  addressed  in  this  study.  An  airborne 
aircraft  will  be  exposed  to  this  radiation  at  cloud  stabili¬ 
zation  time.  As  the  aircraft  flies  toward  the  cloud  the 
radiation  will  increase  to  a  maximum  at  the  horizontal  center 
of  the  cloud  and  will  then  slowly  decrease  as  the  aircraft 
flies  away  from  the  cloud.  There  are  two  major  assumptions 
made  for  this  model: 

1.  The  mean  free  path  of  the  gamma  radiation  is  less 
than  0.1  ox. 

2.  The  aircraft  is  more  than  3  ctx  away  from  ground  zero 
at  t  plus  6  minutes. 

These  assumptions  will  be  discussed  in  greater  detail  later 
in  the  text. 

The  reason  for  the  first  assumption  can  be  explained  as 
follows.  Consider  the  normalized  horizontal  distribution 
function,  f(x,t).  This  function  can  be  rewritten  as 


f(x,t)  = 


-1_  e^( 
}ir  ov  \ 


VAC  » ta-VAC  *  t \ 2 


where  VAC  is  the  velocity  of  the  aircraft.  This  function  is 
now  integrated  with  respect  to  time  from  minus  to  plus  infin¬ 
ity.  After  completing,  this  integration  one  obtains 


r 

f(x,t)  =  -±- 

Y 

ft  Am 

£ 

VAC 

m 
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If  the  mfp  of  the  one  MeV  gamma- ray  is  greater  than  0.1  ux 
the  f(x, t)  function  could  not  be  modeled  as  a  true  normal 
distribution  function.  It  follows  that  as  the  gamma  mfp 
increases  the  f(x,t)  curve  would  vary  and  would  then  be  in¬ 
creasingly  divergent  from  a  true  normal  distribution.  This 
model  would  then  be  inaccurate  and  Eq  (24)  could  not  be  used. 
By  only  using  this  model  for  a  gamma  mfp  of  les  than  0.1  ox 
the  inaccuracy  associated  with  this  model  can  be  greatly 
reduced.  The  author  calculated  only  a  1  to  2%  variation  in 
activity  per  meter  when  the  mfp  was  kept  below  this  limit. 

It  can  also  be  noted  that  the  only  altitudes  affected,  for  a 
one  MeV  gamma-ray,  by  the  above  condition  are  those  greater 
than  11,000  meters.  The  computer  code  outputs  an  error  mes¬ 
sage  when  the  mfp  is  more  than  0.1  ox. 

The  second  assumption  allows  the  f(x,t)  function  to  be 
integrated  in  time  from  plus  to  minus  infinity  with  very 
little  error.  The  true  integration  time  would  be  from  approx¬ 
imately  t+6  minutes  to  ta  plus  mission  time  remaining.  The 
minimum  arrival  time  used  in  this  study  was  30  minutes.  This 
time  would  be  equal  to  approximately  80  ox.  The  integration 
times  used  can  be  compared  to  the  true  integration  times  to 
determine  the  error  associated  with  the  integration  used. 
Consider  the  integral 


f  (x,  t') 


+  6 
=  / 
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/2tt  ax 
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where  ta  -  1/2  hour  to  8  hours.  The  maximum  value  Eq  (25) 
has  would  be  at  ta  =  1/2  hour  and  this  is  almost  infinitesi¬ 
mal.  The  contribution  from  t-  =  ta+8  hours  to  plus  infinity 
is  also  proportionally  small.  Figure  4  shows  that  when  the 
aircraft  is  more  than  a  few  gamma  mfp's  from  the  cloud  the 
exponential  attenuation  of  the  radiation,  by  the  air,  would 
significantly  reduce  the  radiation  received.  Spherical  di¬ 
vergence  also  decreases  the  radiation  received  by  the  air¬ 
craft.  Only  when  the  aircraft  is  within  a  few  gamma  mfp's 
is  the  radiation  significant. 

The  sky-shine  dose  rate  is  calculated  by  the  use  of 
spherical  integration,  namely 

D  =  A' '  '  (x,y,  z,  t)/  SI  ill® - sin<t>ded4>ds  (26) 

0  0  0  P  4tts^ 

Ua 

where  — —  is  the  tissue  absorption  coefficient,  is  the 
P 

attenuation  coefficient  of  air  and  s  is  the  radius  of  the 
cloud.  Equation  (26)  assumes  that  A' ' ' (x,y, z, t)  does  not 
vary  as  the  differential  volume  elements  are  integrated  away 
from  the  aircraft's  location.  After  integrating  Eq  (26)  one 
obtains 

6  =  1.6  x  10_11-A"  (y,z,t)  (l-e~Pts)  (27) 

where  VAC  is  the  aircraft's  velocity  and  1.6  x  lO-^1  is  a 
conversion  factor  to  rad-tissue/hr.  By  letting  s  approach 
infinity  the  exponential  term  in  Eq  (27)  approaches  1.  Eq 
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Figure  4.  Cloud  Model 
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(27)  is  the  unit  time  reference  dose  rate  and  must  be  con¬ 
verted  to  real  time  by  the  Way-Wigner  decay  law. 

A  shielding  model  is  developed  in  order  to  attenuate 
the  sky-shine  dose  which  is  transmitted  through  the  aircraft's 
skin.  The  cabin  area  was  modeled  as  a  right  circular  cylinder 
The  ends  of  the  cylinder  were  then  considered  to  be  infinite 
shields.  This  was  done  to  take  into  account  the  additional 
shielding  supplied  by  the  equipment  and  frame  surrounding 
the  cabin  area.  The  percent  transmitted  through  the  lateral 
area  of  the  cylinder  is  then  approximately  92  percent.  The 
1  MeV  gamma-rays  were  further  attenuated  by  the  factor  e-lJtX 
where  is  the  attenuation  coefficient  for  aluminum  (Ref  12) 
and  x  is  the  skin  thickness.  This  reduced  the  transmitted 
gamma-ray  energy  by  an  additional  2  percent.  Therefore,  an 
attenuation  factor  of  0.9  was  taken  times  the  sky-shine  dose. 
Sky-shine  doses  are  shown  in  Tables  II  through  X.  It  can  be 
seen  from  these  Tables  that  the  dose  is  much  more  significant 
at  early  times  and  higher  altitudes.  This  is  due  to  the  fall 
of  the  particles  and  the  decreased  attenuation  of  the  gamma 
radiation  at  higher  altitudes.  Also  note  from  Eq  (27)  that 
as  the  aircraft's  velocity  increases  the  dose  would  decrease. 

Cabin  Ingestion 

The  aircraft  is  assumed  to  fly  completely  through  the 
cloud  in  the  x  direction  sweeping  out  the  activity  at  a  con¬ 
stant  altitude.  There  is  no  variation  in  the  y  direction  as 
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the  aircraft  flies  along  the  y  centerline.  Therefore,  since 
all  of  the  activity  is  swept  out  in  the  x  direction  the  func¬ 
tion  f (x, t) ,  which  is  a  normal  function,  can  also  be  inte¬ 
grated,  as  was  done  for  the  sky-shine  dose,  from  minus  to 
plus  infinity.  Therefore,  f(x,t)  =  1/VAC  for  the  cabin  in¬ 
gestion  dose  computations  too.  Eq  (3)  then  reduces  to 


(y0- W  ■  '■ 


Ar (za, r, ta) dr£_i 


(28) 


This  is  illustrated  in  Figure  5. 

The  cabin  air  enters  through  the  fifth  stage  of  the 
engine  air  compressor  (Ref  19)  producing  a  constant  mass 
flow  rate  into  the  cabin  of  88.2  kg/min,  as  noted  previously. 
The  mass  flow  rate  is  constant  at  all  aircraft  velocities 
and  altitudes  used  in  this  analysis  (Ref  14) .  The  inlet 
area  is 

Inlet  Area  =  v  c  ,W  (29) 

v  ^  pair 


where  Pair  is  the  air  density  at  the  aircraft's  altitude. 

The  total  activity  ingested  into  the  cabin,  as  the  aircraft 
traverses  the  cloud,  is  Eq  (28)  times  Eq  (29) .  The  nominal 
gamma-ray  energy  is  assumed  to  be  1  MeV  (Ref  6) .  Note  from 
Eq  (19)  that  the  activity  will  decrease  as  velocity  increases. 

The  flow  through  the  cabin  could  not  be  exactly  modeled 
due  to  the  complex  behavior  of  the  particles  after  they  enter 
the  cabin.  The  mean  residence  time  of  a  particle  was  computed 
to  be  approximately  5  minutes.  This  would  be  true  only  if 
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the  particles  flowed  smoothly  through  the  cabin.  This  is 
not  true  because  of  the  complexity  of  the  cabin's  interior. 

It  is  assumed  that  all  of  the  particles  are  equally  distri¬ 
buted  in  the  cabin  from  aircraft  arrival  to  mission  comple¬ 
tion  time.  The  cabin  activity  density  would  then  be  the 
total  activity  in  the  cabin  divided  by  the  cabin's  volume. 

The  feasibility  of  filtering  the  radioactive,  dust 
particles  was  examined.  This  was  determined  to  be  infeasible 
due  to  the  size  variation  of  the  particles  involved  (Ref  8) . 
These  particles  are  from  1/2  to  1800  microns  in  radius. 
Filtering  of  these  par.  cles  would  be  unlikely,  since  many 
filters  of  various  sizes  would  need  to  be  placed  in  tandem  to 
filter  out  even  a  fraction  of  the  ingested  particles.  These 
filters  would  need  to  be  changed  almost  continuously  when 
flying  through  the  cloud  (Ref  11) .  Also,  damage  to  the  air¬ 
craft,  caused  by  the  larger  particles,  may  be  a  greater  haz¬ 
ard  than  the  radiation  threat.  This  would  be  especially  true 
at  lower  altitudes.  This  suggests  that  flying  at  lower  alti¬ 
tudes,  in  order  to  reduce  crew  dose,  may  not  be  the  best 
action  taken. 

To  obtain  cabin  dose  rates,  the  cabin's  volume  is  mod¬ 
eled  as  a  right  circular  cylinder.  It  is  assumed  that  self¬ 
shielding  would  be  zero  and  that  none  of  the  radioactive 
particles,  which  entered  the  cabin,  would  be  exhausted  out 
of  the  cabin  during  the  remainder  of  the  aircraft's  mission. 
Neither  of  these  assumptions  are  true.  It  would  take  a  very 


detailed  analysis  of  the  cabin  interior  to  obtain  correct 
shielding  and  flow  rate  parameters.  The  dose  rate  at  the 
center  of  the  cabin,  due  to  the  suspended  dust  particles, 
is  then 


16.027  2.212  2ir  Ua  (r2+z2)  *5 

2  /  S  !  A-lia-e  utlr  1  rdedrdz 

0  OOP  4„(r2+,2) 


where  r  is  the  radius  of  the  cylinder,  z  is  one-half  the 
height  and  A  is  the  activity  of  the  ingested  particles  at 
one  hour.  Due  to  the  complexity  of  Eq  (30)  a  numerical  in¬ 
tegration,  using  Simpsons  approximation,  was  performed.  The 
equation,  after  the  numerical  integration,  is 

D  =  3.00  Aiia  (31) 

P 

Figure  6  shows  the  cabin  activity  for  a  1  Mt  burst  at 
altitudes  ranging  from  500  to  12,000  meters  and  at  arrival 
times  of  1/2  to  8  hours.  Note  that  the  activity  decreases 
as  the  aircraft's  altitude  decreases  and  for  later  arrival 
times.  This  suggests  that  the  aircraft's  altitude  could.be 
decreased  to  decrease  the  dose  received  by  an  aircrew  member. 
Tables  II  through  X  show  the  cabin  dose  results  from  a  1  Mt 
burst  at  varying  arrival  times  and  altitudes. 
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TABLE  II 


Single  Burst  Dose  Summary  for 
a  lMt  Burst,  TA=  0.5  Hrs,  TR=  8.0  Hrs 


Altitude 

(meters) 

Cabin  Dose 

(rem) 

Sky-Shine 

Dose 

(rem) 

Total 

Dose 

(rem) 

Prominent 

Particle 

(microns) 

12.000 

10.455 

23.737 

34.192 

52.7 

11,500 

9.071 

20.595 

29.666 

62.5 

11,000 

7.690 

17.458 

25.148 

71.3 

10, 500 

6.474 

14.698 

21.172 

81.4 

10,000 

5.364 

12.178 

17.542 

93.5 

9500 

4.394 

9.976 

14.370 

103.0 

9000 

3.578 

8.124 

11.702 

113.0 

8500 

2.910 

6.607 

9.517 

126.0 

8000 

2.378 

5.400 

7.778 

132.0 

7500 

1.958 

4.444 

6.402 

148.0 

7000 

1.627 

3.694 

5.321 

156.0 

6500 

1.368 

3.106 

4.474 

166.0 

6000 

1.141 

2.590 

3.731 

176.0 

5500 

0.982 

2.230 

3.212 

188.0 

5000 

0.857 

1.546 

2.803 

201.0 

4500 

0.753 

1.709 

2.462 

215.0 

4000 

0.666 

1.517 

2.183 

231.0 

3500 

0.591 

1.342 

1.933 

249.0 

3000 

0.526 

1.195 

1.721 

270.0 

2500 

0.468 

1.062 

1.530 

270.0 

2000 

0.416 

0.943 

1.359 

294.0 

1500 

0.368 

0.837 

1.205 

323.0 

1000 

0.325 

0.737 

1.062 

323.0 

500 

0.285 

0.648 

0.933 

357.0 
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TABLE  III 


Altitude 

(meters) 

Cabin  Dose 

(rem) 

Sky-Shine 

Dose 

(rem) 

Total 

Dose 

(rem) 

Prominent 

Particle 

(microns) 

12,000 

5.062 

6.990 

12.052 

34.7 

11,500 

4.427 

6.113 

10.540 

41.0 

11,000 

3.789 

5.233 

9.022 

46.5 

10,500 

3.226 

4.455 

7.681 

52.7 

10,000 

2.710 

3.742 

6.452 

57.3 

9500 

2.257 

3.117 

5.374 

62.5 

9000 

1.874 

2.587 

4.461 

68.2 

8500 

1.558 

2.152 

3.710 

74.5 

8000 

1.303 

1.799 

3.102 

77.7 

7500 

1.100 

1.519 

2.619 

85.2 

7000 

0.938 

1.295 

2.233 

89.2 

6500 

0.808 

•  1,116 

1.924 

98.0 

6000 

0.693 

0.958 

1.651 

103.0 

5500 

0.609 

0.842 

1.451 

108.0 

5000 

0.540 

9.745 

1.285 

113.0 

4500 

0.481 

0.665 

1.146 

119.0 

4000 

0.431 

0.595 

1.026 

126.0 

3500 

0.387 

0.534 

0.921 

132.0 

3000 

0.348 

0.480 

0.828 

140.0 

2500 

0.314 

0.434 

0.748 

148.0 

2000 

0.284 

0.392 

0.676 

156.0 

1500 

0.257 

0.356 

0.613 

166.0 

1000 

0.234 

0.323 

0.557 

166.0 

500 

0.213 

0.294 

0.507 

176.0 
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TABLE  IV 


Single  Burst  Dose  Summary  for 
a  lMt  Burst,  TA=  2.0  Hrs,  TR=  8.0  Hrs 


Altitude 

(meters) 

Cabin  Dose 

(rem) 

Sky-Shine 

Dose 

(rem) 

Total 

Dose 

(rem) 

Prominent 

Particle 

(microns) 

12,000 

1.966 

1.754 

3.720 

24.6 

11,500 

1.731 

1.544 

3.275 

28.0 

11,000 

1.491 

1.330 

2.821 

31.9 

10,500 

1.278 

1.140 

2.418 

34.7 

10,000 

1.082 

1.082 

2.164 

37.7 

9500 

0.908 

0.810 

1.718 

41.0 

9000 

0.760 

0.678 

1.438 

44.6 

8500 

0.637 

0.568 

1.205 

46.5 

8000 

0.538 

0.480 

1.018 

50.5 

7500 

0.458 

0.409 

0.867 

52.7 

7000 

0.395 

0.352 

0.747 

55.0 

6500 

0.345 

0.308 

0.653 

57.3 

6000 

0.301 

0.268 

0.569 

62.5 

5500 

0.268 

0.239 

0.507 

65.3 

5000 

0.241 

0.215 

0.456 

68.2 

4500 

0.219 

0.195 

0.414 

71.3 

4000 

0.199 

0.177 

0.376 

74.5 

3500 

0.182 

0.162 

0.344 

77.7 

3000 

0.166 

0.148 

0.314 

81.4 

2500 

0.152 

0.136 

0.288 

85.2 

2000 

0.140 

0.125 

0.265 

89.2 

1500 

0.128 

0.115 

0.243 

93.5 

1000 

0.118 

0.105 

0.223 

93.5 

500 

0.109 

0.097 

0.206 

98.0 
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TABLE  V 


Single  Burst  Dose  Summary  lor 
a  lMt  Burst,  TA=  3.0  Hrs,  TR=  8.0  Hrs 


Altitude 

(meters) 

Cabin  Dose 

(rem) 

Sky-Shine 

Dose 

(rem) 

Total 

Dose 

(rem) 

Prominent 

Particle 

(microns) 

12,000 

1.019 

0.729 

1.748 

19.5 

11,500 

0.900 

0.644 

1.544 

22.5 

11,000 

0.779 

0.557 

1.336 

24.6 

10,500 

0.670 

0.480 

1.150 

28.0 

10,000 

0.570 

0.407 

0.977 

30.5 

9500 

0.480 

0.343 

0.823 

31.9 

9000 

0.404 

0.289 

0.693 

34.7 

8500 

0.340 

0.243 

0.583 

36  „  2 

8000 

0.288 

0.206 

0.494 

39.3 

7500 

0.246 

0.176 

0.422 

41.0 

7000 

0.212 

0.152 

0.364 

42.8 

6500 

0.185 

0.132 

0.317 

44.6 

6000 

0.161 

0.115 

0.276 

48.4 

5500 

0.144 

0.103 

0.247 

50.5 

5000 

0.129 

0.093 

0.222 

52.7 

4500 

0.117 

0.084 

0.201 

55.0 

4000 

0.107 

0.077 

0.184 

55.0 

3500 

0.098 

0.072 

0.170 

57.3 

3000 

0.090 

0.065 

0.155 

59.8 

2500 

0.083 

0.060 

0.143 

62.5 

2000 

0.077 

0.055 

0.132 

65.3 

1500 

0.072 

0.051 

0.123 

68.2 

1000 

0.066 

0.047 

0.113 

71.3 

500 

0.062 

0.044 

0.106 

74.5 
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TABLE  VI 


Single  Burst  Dose  Summary  for 
a  lMt  Burst,  TA=  4.0  Hrs,  TR=  8.0  Hrs 


Altitude 

(meters) 

Cabin  Dose 

(rem) xl0"l 

Sky-Shine 

Dose 

(rem) xlO-^ 

Total 

Dose 

(rem) 

Prominent 

Particle 

(microns) 

12,000 

6.332 

39.407 

1.027 

17.0 

11,500 

5.609 

34.904 

0.910 

19.5 

11,000 

4.863 

30.262 

0.789 

21.5 

10,500 

4.200 

26.136 

0.681 

23.5 

10,000 

3.580 

22.279 

0.581 

25.7 

9500 

3.027 

18.836 

0.491 

28.0 

9000 

2.551 

15.876 

0.414 

29.3 

8500 

2.153 

13.398 

0.349 

30.5 

8000 

1.828 

11.376 

0.297 

33.2 

7500 

1.565 

9.743 

0.254 

34.7 

7000 

1.354 

8.428 

0.220 

36.2 

6500 

1.183 

7.363 

0.192 

37.7 

6000 

1.034 

6.437 

0.168 

39.3 

5500 

0.921 

5.729 

0.149 

41.0 

5000 

0.827 

5.147 

0.134 

42.8 

4500 

0.749 

4.663 

0.121 

44.6 

4000 

0.682 

4.245 

0.111 

46.5 

3500 

0.623 

3.880 

0.101 

48.4 

3000 

0.573 

3.565 

0.096 

50.5 

2500 

0.528 

3.287 

0.086 

52.7 

2000 

0.489 

3.041 

0.079 

52.7 

1500 

0.454 

2.823 

0.074 

55.0 

1000 

0.423 

2.630 

0.069 

57.3 

500 

0.394 

2.450 

0.064 

59.8 
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TABLE  VII 


Single  Burst  Dose  Summary  for 
a  lMt  Burst,  TA=  5.0  Hrs,  TR=  8.0  Hrs 


Altitude 

(meters) 

Cabin  Dose 

(rem) xl0“^ 

Sky-  Shine 
Dose 

(rem) xlO-2 

Total 

Dose 

(rem) 

Prominent 

Particle 

(microns) 

12,000 

42.252 

23.854 

0.661 

15.4 

11,500 

37.509 

21.176 

0.587 

17.0 

11,000 

32.603 

18.407 

0.511 

19.5 

10,500 

28.221 

15.933 

0.441 

21.5 

10,000 

24.110 

13.611 

0.377 

22.5 

9500 

20.419 

11.528 

0.319 

24.6 

9000 

17.245 

9.736 

0.270 

25.7 

8500 

14.586 

8.235 

0.228 

28.0 

8000 

12.406 

7.004 

0.194 

29.3 

7500 

10.635 

6.004 

0.166 

30.5 

7000 

9.215 

5.203 

0.144 

31  9 

6500 

8.067 

4.555 

0.126 

33.2 

6000 

7.059 

3.985 

0.111 

34.7 

5500 

6.292 

3.552 

0.098 

36.2 

5000 

5.667 

3.200 

0.089 

37.7 

4500 

5.125 

2.894 

0.080 

39.3 

4000 

4.668 

2.635 

0.073 

41.0 

3500 

4.262 

2.406 

0.067 

42.8 

3000 

3.908 

2.206 

0.061 

44.6 

2500 

3.595 

2.030 

0.056 

44.6 

2000 

3 .  314 

1.871 

0.052 

46.5 

1500 

3.0/1 

1.734 

0.048 

48.4 

1000 

2.851 

1.612 

TT 

O 

• 

o 

50.5 

500 

2.656 

1.499 

0.042 

50.5 
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TABLE  VIII 

Single  Burst  Dose  Summary  for 
a  lMt  Burst,  TA=  6.0  Hrs,  TR=  8.0  Hrs 


Altitude 

(meters) 

Cabin  Dose 

(rem)xlO-3- 

Sky-Shine 

Dose 

(rem) xlO"2 

Total 

Dose 

(rem) 

Prominent 

Particle 

(microns) 

12,000 

29.773 

15.631 

0.454 

14.0 

11,500 

26.483 

13.904 

0.404 

15.4 

11,000 

23.066 

12.110 

0.352 

17.9 

10,500 

20.000 

10.500 

0.305 

19.5 

10,000 

17.119 

8.988 

0.261 

20.5 

9500 

14.531 

7.629 

0.222 

22.5 

9000 

12.290 

6.453 

0.187 

23.5 

8500 

10.404 

5.462 

0.159 

24.6 

8000 

8.858 

4.650 

0.135 

26.8 

7500 

7.604 

3.992 

0.116 

28.0 

7000 

6.589 

3.459 

0.101 

29.3 

6500 

5.774 

3.031 

0.088 

30.5 

6000 

5,061 

2.657 

0.077 

31.9 

5500 

4.514 

2.370 

0.069 

33.2 

5000 

4.065 

2.134 

0.062 

34.7 

4500 

3.683 

1.934 

0.056 

34.7 

4000 

3.357 

1.763 

0.051 

36.2 

3500 

3.067 

1.610 

0.047 

37.7 

3000 

2.810 

1.475 

0.043 

39.3 

2500 

2.586 

1.358 

0.039 

41.0 

2000 

2.381 

1.250 

0.036 

42.8 

1500 

2.203 

1.157 

0.034 

42.8 

1000 

2.043 

1.073 

0.031 

44.6 

500 

1.898 

0.997 

0.029 

46.5 
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TABLE  IX 


Single  Burst  Dose  Summary  for 
a  lMt  Burst,  TA=  7.0  Hrs,  TR=  8.0  Hrs 


Altitude 

(meters ) 

Cabin  Dose 

(rem)  xlO"V 

Sky-Shine 

Dose 

(rem) xlO”2 

Total 

Dose 

(rem) 

Prominent 

Particle 

(microns) 

12, 000 

21.872 

10.852 

0.327 

12.6 

11, 500 

19.491 

9.671 

0.292 

14.7 

11,000 

17.003 

8.436 

0.254 

16.2 

10, 500 

14.776 

7.331 

0.221 

17.9 

10,000 

12.670 

6.286 

0.190 

19.5 

9500 

10.772 

5.343 

0.161 

20.5 

9000 

9.120 

4.525 

0.137 

21.5 

8500 

7.732 

3.837 

0.116 

23.5 

8000 

6.588 

3.269 

0.099 

24.6 

7500 

5.655 

2.806 

0.084 

25.7 

7000 

4.905 

2.434 

0.072 

26.8 

6500 

4.298 

2.133 

0.064 

28.0 

6000 

3.766 

1.869 

0.056 

29.3 

5500 

3.362 

1.668 

0.050 

30.5 

5000 

3.028 

1.503 

0.045 

31.9 

4500 

2.749 

1.364 

0.041 

31.9 

4000 

2.505 

1.243 

0.037 

33.2 

3500 

2.291 

1.137 

0.034 

34.7 

3000 

2.103 

1.043 

0.031 

36.2 

2500 

1.935 

0.960 

0.029 

37.7 

2000 

3.784 

0.885 

0.027 

37.7 

1500 

1.648 

0.818 

0.025 

39.3 

1000 

1.528 

0.758 

0.023 

41.0 

500 

1.419 

0.704 

0.021 

41.0 
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TABLE  X 


a  lMt  Burst.,  TA=  8.0  Hrs,  TR=  8.0  Hrs 


Altitude 

(meters ) 

Cabin  Dose 

(rem)  xl0”l 

Sky-Shine 

Dose 

(rem) xlO-2 

Total 

Dose 

(rem) 

Prominent 

Particle 

(microns) 

12, 000 

16.608 

7.875 

0.245 

11.9 

11,500 

14.825 

7.029 

Q .  218 

13.3 

11,000 

12.954 

6.142 

0.191 

15.4 

10, 500 

11.277 

5.347 

0.166 

16.2 

10,000 

9.687 

4.593 

0.143 

17.9 

9500 

8.249 

3.911 

0.122 

18.7 

9000 

6.993 

3.316 

0.103 

20.5 

8500 

5.936 

2.814 

0.087 

21.5 

8000 

5.061 

2.400 

0.075 

22.5 

7500 

4.346 

2.060 

0.064 

23.5 

'’000 

3.770 

1.787 

0.055 

24.6 

o500 

3.303 

1.566 

0.049 

25.7 

6000 

2.895 

1.373 

0.043 

26.8 

5500 

2.583 

1.225 

0.038 

28.0 

5000 

2.326 

1.103 

0.034 

29.3 

4500 

2.111 

1.001 

0.031 

30.5 

4000 

1.924 

0.912 

0.028 

30.5 

?c00 

1.761 

0.835 

0.026 

31.9 

3000 

1.620 

0.768 

0.024 

33.2 

2500 

1.490 

0.707 

0.022 

34.7 

2000 

1.377 

0.653 

0.020 

34.7 

1500 

1.272 

0.603 

0.019 

36.2 

1000 

1.180 

0.559 

0.017 

37.7 

500 

1.094 

0.519 

0.016 

37.7 

IV.  Multiple  Bursts 


The  computer  code  was  initially  developed  to  handle 
only  a  single  burst.  It  was  then  extended  to  compute  doses 
caused  by  multiple  bursts.  The  methods  of  Hickman  (Ref  10) 
and  Crandley  (Ref  4)  were  used  to  develop  this  extension. 

A  field  of  150  km  by  150  km  with  300  bursts  spaced  pro¬ 
portionally  on  the  field  is  used  in  this  study.  These  values 
can  easily  be  changed.  The  individual  clouds  are  assumed  to 
overlap  each  other.  This  is  true  when  ox>L/nx  and  Oy>W/ny. 
Where  nx«ny=N,  the  number  of  bursts,  L=field  length  and 
W=field  width.  Figure  7  illustrates  this  concept.  A  calcu¬ 
lation  is  made  at  arrival  time  to  insure  that  ox>L/nx  and 
cry>W/ny.  The  activity  density  of  the  cloud  is 


A' ' ' (x,y, z, ta)  =  530  X  3.7  x  1016  •  YLD  •  FF  * 
f(x,  ta)  •  f(y,ta)  •  f(z,ta)  fdpsl 

l  m3J 


(32) 


Crandley  calculated  f(y,t)  as 


f(y,ta)  =  ms  r 

u  i  r  a 


150  -75  ^2tt  ay(ta) 


(33) 


which  is  the  sum  of  the  overlapped  Gaussian  distributions. 
The  f(x,  t)  function  is  found  by  the  same  method. 

Bridgman  and  Hickman  integrated  Eq  (33) ,  along  with  its 
x  direction  counterpart,  and  obtained  the  activity  density 
at  the  Cloud  center  at  time  of  arrival  as 
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A"Mx0.y0.*a.ta,  [fi]  (34) 

They  also  found  the  activity  per  area  to  be 

‘"(»o.*.'V  •  i  A;i*a,ta)  [£i]  os) 

This  provides  a  burst  amplification  factor  for  cabin  dose  of 
N/L»2iray  (tg)  .  This  result  is  used  to  compute  the  multiburst 
cabin  dose  in  this  study. 

An  amplification  factor  for  sky-shine  dose  is  found  in 
a  similar  manner.  The  activity  density  for  a  single  burst 
is  equated  to  the  activity  density,  Eq(33),  for  a  multiple 
burst.  A  burst  amplification  factor  of  N/WL* f (x, ta) • f (y , ta) 
was  then  calculated.  The  sky-shine  dose  is  then  computed  for 
a  worst  case  dose  by  integrating  for  the  time  it  took  the 
aircraft  to  fly  over  3  times  the  standard  deviation  of  the 
cloud.  The  aircraft  would  experience  the  activity  density  at 
cloud  center  during  this  time. 

The  results  for  the  multiple  burst  case  are  summarized 
in  Tables  XI  through  XIX.  Figure  8  is  the  cabin  activity 
versus  altitude,  for  arrival  times  of  1/2  to  8  hours,  for  the 
multiple  burst  case. 
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/*.»  TABLE  XI 

& 

Multiple  Burst  Dose  Summary  for 
300  One-Mt  Bursts.  TA=  0.5  Hrs,  TR=  8.0  Hrs 


Altitude 

(meters) 

Cabin  Dose 

(rem) 

Sky-Shine 

Dose 

(rem) 

Total 

Dose 

(rem) 

Prominent 

Particle 

(microns) 

12,000 

301.969 

677.4/7 

979.446 

52.7 

11, 500 

261.995 

587.795 

349.790 

62.5 

11,000 

222.084 

498.254 

720.338 

71.3 

10,500 

186.977 

419.489 

606.466 

81.4 

10,000 

154.922 

347.573 

502.495 

93.5 

9500 

126.904 

284.713 

411.617 

J03.0 

9000 

103.351 

231.872 

335.223 

113.0 

8500 

84.054 

188.578 

272.632 

126.0 

8000 

68.689 

154.106 

222.800 

137.0 

7500 

56.537 

128.843 

185.380 

148.0 

7000 

46.997 

105.438 

152.435 

156.0 

6500 

, 39.514 

88.650 

128.164 

166.0 

6000 

32.943 

73.910 

106.853 

176.0 

5500 

28.374 

63.658 

92.032 

188.0 

5000 

24.750 

55.528 

80.278 

201.0 

4500 

21.747 

48.790 

20.537 

215.0 

4000 

19.233 

43.149 

62.382 

231.0 

3500 

17.068 

38.292 

55.360 

249.0 

3000 

15.196 

34.094 

49.290 

270.0 

2500 

13.513 

30.316 

43.829 

270.0 

2000 

12.003 

26.928 

38.931 

294.0 

1500 

10.644 

23.879 

34.523 

323.0 

1000 

9.382 

21.049 

30.431 

323.0 

500 

8.243 

18.494 

26.737 

357.0 
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TABLE  XII 


Multiple  Burst  Dose  Summary  for 
300  One-Mt  Bursts.  TA=  1.0  Hr.  TR=  8.0  Hrs 


Altitude 

(meters ) 

Cabin  Dose 

(rem) 

Sky-Shine 

Dose 

(rem) 

Total 

Dose 

(rem) 

Prominent 

Particle 

(microns) 

12,000 

186.074 

272.788 

458.862 

34.7 

11,500 

162.732 

238.568 

401.300 

41.0 

11,000 

139.292 

204.205 

343.497 

46.5 

10,500 

118.594 

173.861 

292.455 

52.7 

10,000 

99.620 

146.046 

245.666 

57.3 

9500 

82.968 

121.633 

204.601 

62.5 

9000 

68.875 

100.973 

169.848 

68.2 

8500 

57.275 

83.967 

141.242 

74.5 

8000 

47.893 

70.212 

118.105 

77.7 

7500 

40.426 

59.266 

99.692 

85.2 

7000 

34.481 

50.550 

85.031 

89.2 

6500 

29.707 

43.551 

73.258 

98.0 

6000 

25.493 

37.374 

62.867 

103.0 

5500 

22.403 

32.843 

55.246 

108.0 

5000 

19.841 

29.088 

48.929 

113.0 

4500 

17.696 

25.943 

43.639 

119.0 

4000 

15.834 

23.213 

39.047 

126.0 

3500 

14.220 

20.846 

35.066 

132.0 

3000 

12.788 

18.748 

31.536 

140.0 

2500 

11.549 

16.931 

28.480 

148.0 

2000 

10.439 

15.304 

25.743 

156.0 

1500 

9.465 

13.875 

23.340 

166.0 

1000 

8.588 

12.590 

21.178 

166.0 

500 

7.833 

11.483 

19.316 

176.0 
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TABLE  XIII 


Multiple  Burst  Dose  Summary  for 
300  One-Mt  Bursts,  TA=  2.0  Hrs,  TR=  8.0  Hrs 


Altitude 

(meters) 

Cabin  Dose 

(rein) 

Sky-Shine 

Dose 

(rem) 

Total 

Dose 

(rem) 

Prominent 

Particle 

(microns) 

12/ 000 

105.868 

106.006 

211.874 

24.6 

11/500 

93.170 

93.291 

186.461 

28.0 

11,000 

80.252 

80.357 

160.609 

31.9 

10,500 

68.817 

68.907 

137.724 

34.7 

10,000 

58.231 

58.306 

116.537 

37.7 

9500 

48.867 

48.932 

97.799 

41.0 

9000 

40.892 

40.945 

81.837 

44.6 

8500 

34.292 

34.337 

68.629 

46.5 

8000 

28.947 

28.987 

57.934 

50.5 

7500 

24.663 

24.695 

49.358 

52.7 

7000 

21.269 

21,297 

42.566 

55.0 

6500 

18.565 

18.589 

37.154 

57.3 

6000 

16.189 

16.210 

32.399 

62.5 

5500 

14.421 

14.440 

28.861 

65.3 

5000 

12.997 

13.014 

26.011 

68.2 

4500 

11.778 

11.794 

23.572 

71.3 

4000 

10.712 

10.726 

21.438 

74.5 

3500 

9.777 

9.790 

19.567 

77.7 

3000 

8.949 

8.961 

17.910 

81.4 

2500 

8.198 

8.209 

16.407 

85.2 

2000 

7.527 

7.536 

15.063 

89.2 

1500 

6.919 

6.928 

13.847 

93.5 

1000 

6.369 

6.377 

12.746 

93.5 

500 

5.867 

5.875 

11.742 

98.0 
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TABLE  XIV 

1 

Multiple 

Burst  Dose  Summary  for 

.Si 

•3 

300  One-Mt  Bursts,  TA=  3.0 

Hrs ,  TR=  8.0 

Hrs 

1 

1  • 

.  j 

j 

■ii 

Altitude 

(meters) 

Cabin  Dose 

(rem) 

Sky-Shine 

Dose 

(rem) 

Total 

Dose 

(rem) 

Prominent 

Particle 

(microns) 

i 

12,000 

73.004 

60.474 

133.478 

19.5 

s 

11, 500 

64.489 

53.420 

117.909 

22.5 

11,000 

55.771 

46.199 

101.970 

24.6 

fj 

10,500 

48.028 

39.784 

87.812 

28.0 

a 

10,000 

40.806 

33.802 

74.608 

30.5 

14 

9500 

34.398 

28.494 

62.892 

31.9 

:j 

•-) 

9000 

28.906 

23.945 

52.851 

34.7 

& 

.j 

8500 

24.3.33 

20.157 

44.490 

36.2 

s  ^ 

8000 

7500 

20.603 

17.599 

17.067 

14.578 

37.670 

32.177 

39.3 

41.0 

J 

.1 

;f . 

7000 

15.183 

12.577 

27.760 

42.8 

\ 

i! 

6500 

13.254 

10.979 

24.233 

44.6 

6000 

11.561 

9.577 

21.138 

48.4 

5500 

10.298 

8.531 

18.829 

50.5 

1 

5000 

9.269 

7.678 

16.947 

52.7 

4500 

8.415 

6.971 

15.386 

55.0 

j 

4000 

7.679 

6.361 

14.040 

55.0 

2 

3500 

7.031 

5.824 

12.855 

57.3 

8 

3000 

6.474 

5.363 

11.837 

59.8 

'3 

* 

2500 

5.983 

4.956 

10.939 

62.5 

l 

2000 

5.532 

4.583 

10.115 

65.3 

! 

1500 

5.128 

4.248 

9.376 

68.2 

“i 

j 

1000 

4.759 

3.943 

8.702 

71.3 

s 

i 

500 

4,422 

3.663 

8.085 

74.5 
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TABLE  XV 


Multiple  Burst  Dose  Summary  for 
300  One-Mt  Bursts,  TA=  4.0  Hrs,  TR=  8.0  Hrs 


Altitude 

(meters) 

Cabin  Dose 

(rem) 

Sky-Shine 

Dose 

(rem) 

Total 

Dose 

(rem) 

Prominent 

Particle 

(microns) 

12,000 

54.933 

39.780 

94.713 

17.0 

11, 500 

48.656 

35.235 

83.891 

19.5 

11, vOO 

42.184 

30.548 

72.732 

21.5 

10,500 

36.432 

26.383 

62.815 

23.5 

10,000 

31.056 

22.489 

53.545 

25.7 

9500 

26.257 

19.015 

45.272 

28.0 

9000 

22.131 

16.027 

38.158 

29.3 

8500 

18.677 

13.525 

32.202 

30.5 

8000 

15.858 

11.484 

27.342 

33.2 

7500 

13.581 

9.835 

23.416 

34.7 

7000 

11.749 

8.508 

20.257 

36.2 

6500 

10.264 

7.433 

17.697 

37.7 

6000 

8.973 

6.498 

15.471 

39.3 

5500 

7.986 

5.783 

13.679 

41.0 

5000 

7.175 

5.196 

12.371 

42.8 

4500 

6.501 

4.708 

11.209 

44.6 

4000 

5.917 

4.285 

10.202 

46.5 

3500 

5.408 

3.916 

9.324 

48.4 

3000 

4.969 

3.598 

8.567 

50.5 

2500 

4.582 

3.318 

7.900 

52.7 

2000 

4.240 

3.070 

7.310 

52.7 

1500 

3.935 

2.850 

6.785 

55.0 

1000 

3.666 

2.654 

6.320 

57.3 

500 

3-415 

2.473 

5.888 

59.8 
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TABLE  XVI 


Multiple  Burst  Dose  Summary  for 


300  One-Mt  Bursts,  TA=  5.Q  Hrs ,  TR=  8 . 0  Hrs 


Altitude 

(meters) 

Cabin  Dose 

(rem) 

Sky-Shine 

Dose 

(rem) 

Total 

Dose 

(rem) 

Prominent 

Particle 

(microns) 

12,  000 

43.513 

28.668 

72.181 

15.4 

11,500 

38.627 

25.449 

64.076 

17.0 

11,000 

33.576 

22.121 

55.697 

19.5 

10/ 500 

29.063 

19.148 

48.211 

21.5 

10,000 

24.829 

16.358 

41.187 

22.5 

9500 

21.028 

13.854 

34.882 

24.6 

9000 

17.759 

11.700 

29.459 

25.7 

8500 

15.021 

9.896 

24.917 

28.0 

8000 

12.776 

8.417 

21.193 

29.3 

7500 

10.952 

7.216 

18.168 

30.5 

7000 

9.490 

6.252 

15.742 

31.9 

6500 

8.309 

5.474 

13.783 

33.2 

6000 

7.269 

4.789 

12.058 

34.7 

5500 

6.480 

4.269 

10.749 

36.2 

5000 

5.837 

3.845 

9.662 

37.7 

4500 

5.278 

3.478 

8.756 

39.3 

4000 

4.807 

3.167 

7.974 

41.0 

3500 

4.389 

2.892 

7.281 

42.8 

3000 

4.025 

2.652 

6.677 

44.6 

2500 

3.703 

2,439 

6.142 

44.6 

2000 

3.413 

2.248 

5.661 

46.5 

1500 

3.163 

2.084 

5.247 

46.4 

1000 

2.940 

1.937 

4.077 

50.5 

500 

2.735 

1.802 

4.537 

50.5 
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TABLE  XVII 

Multiple  Burst  Dose  Summary  for 
300  One-Mt  Bursts,  TA=  6.0  Hrs,  TR=  8.0  Hrs 


Altitude 

(meters) 

Cabin  Dose 

(rem) 

Sky-Shine 

Dose 

(rem) 

Total 

Dose 

(rem) 

Prominent 
Particle 
(microns ) 

12/ 000 

35.690 

21.909 

57.599 

14.0 

11/500 

31.746 

19.488 

51.234 

15.4 

11/000 

27.650 

16.973 

44.623 

17.9 

10,500 

23.974 

14.717 

38.691 

19.5 

10,000 

20.521 

12.597 

33.118 

20.5 

9500 

17.418 

10.692 

28.110 

22.5 

9000 

14.733 

9.044 

23.777 

23.5 

8500 

12.471 

7.656 

20.127 

24.6 

8000 

10.618 

6.518 

17.136 

26.8 

7500 

9.115 

5.595 

14.710  - 

--■23  VO - 

7000 

7.899 

4.849 

12.748 

29.3 

6500 

6.921 

4.249 

11.170 

30,5 

6000 

6 . 066 

3.724 

9.790 

31.9 

5500 

5.411 

3.322 

8.733 

33.2 

5000 

4.873 

2.991 

7.864 

34.7 

4500 

4.415 

2.710 

7.125 

34.7 

4000 

4.025 

2.471 

6.496 

36.2 

3500 

3.676 

2.257 

5.993 

37.7 

3000 

3.368 

2.067 

5.435 

39.3 

2500 

3.100 

1.903 

5.003 

41.0 

2000 

2.854 

1.752 

4.606 

42.3 

1500 

2.641 

1.621 

4.262 

42.8 

1000 

2.449 

1.504 

3.953 

44.6 

500 

2.276 

1.397 

3.673 

46.5 
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TABLE  XVIII 


Multiple  Burst  Dose  Summary  for 
300  One-Mt  Bursts,  TA=  7.0  Hrs.  TR=  8.0  Hrs 


Altitude 

(meters ) 

Cabin  Dose 

(rem) 

Sky-Shine 

Dose 

(rem) 

Total 

Dose 

(rem) 

Prominent 

Particle 

(microns) 

12,000 

30.004 

17.431 

47.435 

12.6 

11,500 

26.738 

15.533 

42.271 

14.7 

11,000 

23.325 

13.550 

36.875 

16.2 

10,500 

20.269 

11.775 

32.044 

17.9 

10,000 

17.381 

10.097 

27.478 

19.5 

9500 

14.777 

8.854 

23.631 

20.5 

9000- 

12.510 

7.268 

19.778 

21.5 

"'1*500 

10.607 

6.162 

16.229 

23.5 

8000 

9.038 

5.250 

14.288 

24.6 

7500 

7.757 

4.506 

12.263 

25.7 

7000 

6.728 

3.909 

10.637 

26.8 

6500 

5.896 

3.425 

9.321 

28.0 

6000 

5.167 

3.002 

8.169 

29.3 

5500 

4.612 

2.679 

7.291 

30.5 

5000 

4.154 

2.413 

6.567 

31.9 

4500 

3.771 

2.190 

5.961 

31.9 

4000 

3.436 

1.996 

5.432 

33.2 

3500 

3.142 

1.825 

4.967 

34.7 

3000 

2.884 

1.676 

4.560 

36.2 

2500 

2.654 

1.542 

4.196 

37.7 

2000 

2.447 

1.422 

3.869 

37.7 

1500 

2.261 

1.313 

3.574 

39.3 

1000 

2.096 

1.217 

3.313 

41.0 

500 

1.946 

1.131 

3.077 

41.0 
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TABLE  XIX 


Multiple  Burst  Dose  Summary  for 
300  One-Mt  Bursts.  TA=  8.0  Hrs,  TR=  8.0  Hrs 


Altitude 

(meters) 

Cabin  Dose 

(rem) 

Sky-Shine 

Dose 

(ram) 

Total 

Dose 

(rem) 

Prominent 

Particle 

(microns) 

12,000 

25.702 

14.284 

39.986 

11.9 

11, 500 

22.943 

12.750 

35.693 

13.3 

11,000 

20.048 

11.142 

31.190 

15.4 

10,500 

17.453 

9.699 

27.152 

16.2 

10,000 

14.992 

8.331 

23.323 

17.9 

9500 

12.766 

7.095 

19.861 

18.7 

9000 

10.823 

6.014 

16.837 

20.5 

8500 

9.186 

5.105 

14.291 

21.5 

8000 

7.833 

4.353 

12.186 

22.5 

7500 

6.725 

3.738 

10.463 

23.5 

7000 

5.834 

3.242 

9.076 

24.6 

6500 

5.112 

2.841 

7.953 

25.7 

6000 

4.481 

2.490 

6.971 

26.8 

5500 

3.997 

2.221 

6.218 

28.0 

5000 

3.600 

2.001 

5.601 

29.3 

4500 

3.267 

1.816 

5.083 

30.5 

4000 

2.978 

1.655 

4.633 

30.5 

3500 

2.725 

1.514 

4.239 

31.9 

3000 

2.507 

1.393 

3.900 

33.2 

2500 

2.306 

1.282 

3.588 

34.7 

2000 

2.131 

1.184 

3.315 

34.7 

1500 

1.968 

1.094 

3.062 

36.2 

1000 

1.825 

1.014 

2.839 

37.7 

500 

1.693 

0.941 

2.634 

37.7 
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V.  Conclusions  and  Recommendations 


Conclusions 

The  results  obtained  in  this  study  show  a  decrease  in 
dose  received  of  20  to  40  percent  for  sky-shine  dose  as 
compared  to  Hickman's  worst  case  analysis  for  a  single  burst. 
For  the  sky-shine  dose,  integrating  the  f(x,t)  function  from 
minus  to  plus  infinity,  instead  of  using  Hickman's  worst  case 
approximation,  reduces  the  dose  received  significantly.  Also, 
the  10%  shielding  approximation  contributes  to  this  dose 
reduction.  The  cabin  ingestion  dose  results  for  a  single 
burst  approximate  Hickman's  results.  The  reason  the  cabin 
ingestion  doses  are  close  is  because  Hickman's  approximation 
swept  out  approximately  96*  of  the  particles?  while  this 
study's  analysis  swept  out  all  of  the  particles.  This  still 
amounts  to  a  significant  dose  for  a  multiple  burst  scenario. 

The  dose  received  is  larger  at  earlier  arrival  times  and 
higher  altitudes.  The  dose  received  by  the  aircrew  is  directly 
dependent  on  arrival  time,  aircraft  altitude  and  aircraft 
velocity.  Longer  mission  times,  after  cloud  penetration, 
and  more  bursts,  in  the  multiple  burst  approach,  would  pro¬ 
portionally  increase  the  dose  each  crew  member  receives . 

The  feasibility  of  filtering  radioactive  particles  before 
they  are  able  to  enter  the  cabin  was  previously  discussed. 

It  was  mentioned  that  the  filtering  of  these  particles  was 
infeasible.  Verbal  information  concerning  partial  filtering 
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was  received  by  the  author  (Ref  19) .  The  data  supporting 
this  has  not  been  received  at  the  time  of  this  analysis. 

A  shielding  factor  was  incorporated  for  the  sky-shine 
dose.  This  factor,  although  not  exact,  is  most  likely  less 
than  the  true  shielding  factor.  This  is  because  only  the 
lateral  surface  area  surrounding  the  crew  was  shielded  by 
the  skin  thickness.  Shielding  of  the  end  areas  was  assumed 
to  be  infinite.  This  is  a  trivial  assumption  since  the  end 
areas  contain  less  than  8  percent  of  the  total  area. 

Attenuation  of  the  radiation  by  air  at  varying  distances 
from  the  cloud  was  also  incorporated  into  this  model.  There¬ 
fore,  this  is  not  a  worst  case  approach  to  the  sky-shine  dose. 
The  assumptions  made  for  this  model  were  detailed  in  the  sky- 
shine  section  of  this  report.  It  was  shown  for  the  arrival 
times  and  altitudes  used  that  the  sky-shine  doses  received 
are  very  accurate. 

Recommendations 

There  are  four  recommendations  to  be  made.  First,  a 
single  westerly  wind  was  used  in  this  analysis.  A  very 
detailed  analysis  incorporating  nominal  wind  patterns  for 
various  latitudes,  longitudes  and  altitudes  for  the  continental 
United  States  could  be  made.  These  changes  could  then  be 
included  in  a  revised  computer  code.  Although,  this  is  a 
major  change  and  a  very  difficult  one  it  would  give  wartime 
planners  a  more  accurate  picture  of  the  radioactive  cloud's 
distribution  as  a  function  of  its  position.  Secondly,  a 


more  detailed  analysis  of  the  cabin's  interior  shielding  and 
flow  rate  parameters  could  be  made  to  predict  a  more  accurate 
cabin  dose.  Due  to  the  complicated  nature  of  the  cabin's 
interior,  an  interior  cabin  simulator  would  probably  be  needed 
to  accomplish  this.  Plow  rate  parameters,  size  and  percent 
of  particles  settling  out  and  rate  and  size  of  particles  ex¬ 
hausted  could  then  be  found  empirically.  Thirdly,  a  detailed 
analysis  of  the  materials  and  their  respective  thicknesses 
surrounding  the  crew  would  need  to  be  done  to  obtain  a  true 
shielding  factor  for  the  sky-shine  radiation.  Fourth,  there 
is  a  data  source  that  could  be  used  to  bypass  the  assumption 
that  all  the  gamma- rays  have  an  energy  of  1  MeV  (Ref  6) . 
Therefore,  depending  on  the  time  elapsed  after  the  burst,  a 
more  accurate  gamma-ray  energy  could  be  used.  This  would 
cause  the  parameters  which  are  dependent  on  gamma-ray  energy, 
such  as  tissue  absorption  and  attenuation  coefficients,  to 
change  with  time. 
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Appendix  A 


PF  =  Fission  fraction 
YLD  *  Yield  in  kilotons 
ACR  *  Aircraft  range  to  cloud 
VAC  a  True  airspeed 
ACA  =  Aircraft  altitude  in  meters 
TR  s  Mission  time  remaining 
R  *  Particle  radii 
TA  =  Arrival  time 
VOL  a  cabin  volume 
HR  *  Arrival  time  in  hours 
W  a  Mass  flow  rate 
HC  =  Cloud  center  height 

Cl  =  Intercept  for  initial  particle  height  equation 
C2  =  Slope  for  initial  particle  height  equation 
ANAX  =  Cloud  activity 
RP  *  particle  density 
PG  =  Sea  level  air  pressure 
TC  *  Time  constant 
SO  «  Sigma  naught 
TT  =  True  time  "sed 
SX  =  Sigma  x 
FX  =  f  (x,  t) 


DALT  =  Delta  altitude 


DT  *  Delta  time 

ATOT  »  Total  activity  at  a  specific  altitude 
SQ  »  Sigma  z  at  the  aircraft  altitude 
Z  «  Particle  distribution  center  height 
T  =  Temperature 
P  =  Pressure 
RA  =  Air  density 
DV  »  Dynamic  viscosity 
KV  *  Kinematic  viscosity 

RC  =  Reynolds  number  squared  times  the  drag  coefficient 
RE  =  Reynolds  number 
VP  *  Particle  velocity 
T1  =  Time  step 
AR  «  Activity  per  meter 
PER  *  Percent  of  total  activity 
PARMAX  *  Dominent  particle's  activity 
RADMAX  *  Dominent  particle's  radius 
FY  a  f  (y,  t) 

SY  *  sigma  yv 

A2  a  Activity  per  meter  squared 
A3  =  Activity  per  meter  cubed 
DENS  =  Air  density 
CACT  =  Cabin  activity 
DRATE  a  Dose  rate 

TD  *  Cabin  tissue  dose 


6C 


MUA  =  Tissue  absorption  coefficient 
MUT  =  Attenuation  coefficient  in  air 
D1  =  Dose  rate  for  sky-shine  dose 
SD  =  Sky-shine  dose 
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Appendix  B 
Dose  Program 

A  brief  explanation  of  the  dose  program  and  its  listing 
are  contained  in  this  appendix.  This  program  can  be  run 
either  interactively  or  by  card  deck.  Any  initial  parameter 
can  be  changed,  although  there  are  four  parameters  that  are 
most  likely  to  be  changed.  These  parameters  are: 

1.  Yield  in  Kilotons  (YLD) 

2.  Aircraft  range  to  burst  in  meters  (ACR) 

3.  True  airspeed  in  meters  per  second  (VAC) 

4.  Mission  time  remaining  in  hours  (TR) 

Initially,  the  program  calculates  and  displays  the  time 
of  arrival,  cloud  center  height  and  the  aircraft's  distance 
to  the  burst.  It  then  calculates  and  displays,  for  24  alti¬ 
tudes  from  12,000  to  500  meters,  activity  density,  activity 
area,  aircraft  altitude,  cabin  dose,  gamma-ray  mean  free  path 
and  sigma  x.  The  last  two  values  are  compared  with  each 
other  to  insure  an  accurate  calculation  of  the  sky-shine  dose. 
Dominent  particle  size  and  the  percent  of  the  total  activity 
it  contributes  are  also  displayed. 

Another  version  also  calculates  and  displays  the  particle 
distribution  center  height  and  the  activity  per  meter  at  this 
altitude  for  each  particle. 
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MRISRT  PARTICLE  SIZE  IS  .MM743  II  RMIOS 

THIS  PARTICLE  COITIIMTES  J.S4341J74472H  S  TOTAL  ACTIVITY 

FM  PARTICLES  KTSER  4.S-7  MR  .000144  H  MAIN 

ACTIVITT/R2  -  *2.333704*340  CI/R2 

MTIVm/R3  •  .0Y27Y45Y741Y72  Cl/IO 

AIOOFT  STITOE  •  Ml.  STEW 

CMII  ACTIVITY  «  1.711191713217  Cl 

Tim  MS  *  COHOE  ♦  0.  HOI  It  1.930170447343  MI8 

MS  HIE  *  .I74344043IY73  OAOt-TIMS/HOR 

KYMIK  NS  >  2.131904*0413  RMS 

COR  0*  •  317.1943344140  0  1100  I  -  4224.012431774  0 

mttmtmtmttmtmmtmtmttttmmmmtm 


tmmmmtmtmtmmtmtmtmmmttmmm 


otonar  pmticle  its  is  .mmn  o  amirs 


Ivi  n  .  ..  R  ■  i 


0.447  CP  SCOOO  COOtLWNO  THE. 

nn  particle  CMnionis  4.12124403104*  m  nrs  actiyity 

OH  POmOES  STEER  4.S-7  00  .000144  o  OM IS 

ACTIYITY/B  -  171.49079734  CI*2 

ACTIYITT/W3  -  .M23NIA137744  CI/H3 

AIDOH  STITME  -  MM.  STBS 

CMIO  ACTIYITY  •  0. 137104121741  Cl 

Tim  MS  n  EIPIOOE  ♦  0.  MSS  II  1.302412433241  RAM 

ME  AAIE  ■  .732*4942M94  MO-T19K/M 

SnSOIK  MB-  1.7*110441134  RAN 

OHO  VP  ■  2*. 0074074034  fl  1100  I  •  4224.0134317*  It 

tmtmmimtmttimmmmttmitmttttmttm 


mtmtmmmmttmmtittmmtmimtmttttttt 


IMNERT  NOTICLE  HIE  II  .MHB2  0  AMIS 

mi  NOTICLE  COmiMIES  4.447Y727043MX  S  TITS  ACTIVITY 

FM  PAfTICLEt  STSEM  4.S-7  MO  .011174  I  OMIM 


*wm/jc  •  779.419798191  ci/n2 
activity/b  >  .173412131312*4  ci/hj 

AlCAAFT  ALTITME  •  73M.  lETOS 

CMIH  ACTIVITY  -  4.848478121123  Cl 

TISMI  UK  TV  EIP9MK  ♦  I.  MRS  IS  I.899792MII3I  IMS 

USE  MTT  •  .418543728911  MM-TISSK/HM 

SKTaUK  MX  ■  1.31844143374  MM 

MM  KP  •  2I2.I224234IM  8  SIM  I  ■  4224.112431774  8 


tmtmmiimiitsmtttmititttmittmtittmiiitti 


NRIKKT  PMTICLE  SIZE  IS  .NMIV2  8  MIMS 

MIS  PMTICLE  CNTtlNTES  4.97*2113144891  V  TVT4L  ACTlYin 

FOR  PARTICLES  KMEER  4.K-7  AM  .MUM  R  RMIM 

ACTIVITY/R2  •  784.2I843N43  Cl/83 

ACTIVITY/®  •  .#4*3187*497433  CI/R3 

AICRKT  ALTITME  •  7M8.  KTERS 

CMIH  ICTIVin  •  S.BMMM2M3  Cl 

TUSK  NK  Tl  EIPNNK  ♦  L  RHH  II  .93M417127323  IMS 

memie*  .vmmnm  mm-timk/rmr 

KYSRIK  HK  •  1.293314139984  MM 

SMR  m  -  244.399729M94  R  SNM  I  •  4211.11243177*  8 

ttmttmttttimittinmiitMittitmiiiitmmttttt 


ttttitttttmttttmtmimtiitstmtmttttimitttmi 


Minor  particle  size  is  .nmm  r  nmim 

ms  PARTIOE  CMTRIMTES  3.4S4HS241243Z  K  TVTAL  ICTIVin 

FM  PARTICLES  SOKE*  4.K-7  AM  .NUN  I  RMIM 

ACTtVITY/M  -  442.121 124411  CI/R2 

RCT1VITV/R3  •  .14143444324313  CI/K3 

Main  ALTITME  •  43M.  MIMS 

CAM  ACTIVITY  -  3.847191244123  Cl 

TIMK  MX  I*  EVMME  ♦  I.  MMS  IS  .NH44N33334  MM 

NX  MTf  >  .4343289947771  ARN-T1IMI/MM 

KYSRIK  HK  •  1. 1 19938434734  MM 

•MR  8FV  >  231.74528N22I  H  SIMA  I  >  4224.812431774  R 

iiiimtiiitiismtmiitmttttttitisttmtitmimtti 


iitMtititmttiMitmmtitmtmtttttimimimmi 


NRM9T  PARTICLE  SIS  IS  .MM3  R  RMIM 


74 


THIS  PANTICLE  CMWIUTES  4. 911/933423 1 71  OF  TOTAL  ACTIVITY 
t(K  PARTICLES  K TWEEN  .999914  AM  .999291  A  RA9IUS 

ACTIVITY/N2  >  592.4921 127494  C1/K2 
ACTIVITY/H3  *  .93392434441333  CI/M3 
AICRAFT  ALTITOK  *  4999.  UTERS 
CABIN  ACTIVITY  >  4. 33 I 3449331 I 3  C! 

TISSUE  lOSt  TO  UPOSURE  ♦  I.  HOMS  II  .493^33439044  MIS 

MSE  HATE  «  .3999421199934  M9S-TI8M/NM 

siXYSMINE  MSE  «  .9574744394434  RAM 

BANNA  HFP  »  239.9737217274  N  SIMA  I  •  4224.912431774  N 

tMtitiiiiiiiiiiiiimittiiiiiittitittiiittiittiiiiitttt 


tmmttmttimtmttttmttmittmtttttitttttitmi 


•ONTNENT  PARTICLE  SIZE  IS  .999191  R  NA9I9S 

THIS  PARTICLE  CWTRIMTES  4.402199411741  W  TOTAL  RCTIVITV 

FOR  PARTICLES  KTKEN  .9999233  AM  .999213  N  RA9IM 

ACTIVITY /N2  -  349.9192934743  C1/R2 

ACTIVFTY/R3  *  .931994391411  Cl m 

AICMFT  ALTITOE  -  3399.  NETOS 

CABIN  ACTIVITY  -  3. 3B43V99V3ZX7  Cl 

T1SME  BNE  11  OWE  ♦  I.  MMS  II  .40M7IMB23  MM 

ME  NATE  «  .3477M337R979  MM-TOTNE/MM 

BRVBI1M  ME  •  .0413174243333  MR 

9MM  Wf  ■  225.3V77f(f3699  N  SIMA  1  -  4224.912431774  N 

tttmmtmtmmtimimiitiiittmtiimittitttm 


tiiiiiiiiiimimiiiitittmtitmiiitiiitMttifiitmM 


BMMERT  PARTICLE  SIZE  IS  .999113  H  MBIM 

THIS  PARTICLE  GMTRIfiUTES  4  >221113114771  IF  TOTAL  ACTIYITY 

FOR  PARTICLES  BETUEEN  .909.31 I  AM  .990215  I  RM1M 

ACTIVITY/ID  •  393.9157352999  CI/R2 

ACTIVITY/ID  •  .9477IIB2027977  Cl /IT 

AICMFT  ALTITIK  «  3N9.  NETEH 

CAIIH  mYITY  •  3.371944974942  Cl 

TISSUE  C  TO  E1P0MRE  ♦  I.  NHRI  IS  .33977H9497N  MM 

BOSE  BATE  -  .3933*23941041  RMS-TI9SUE/RMR 

SKTSNIM  IOK  >  .7433493993333  MBS 

BANNA  RFP  •  JI3.J74071937B  N  SUM  I  -  4224.912431774  H 

titttiiitimtmtimiiiiiitmtiitifittimimitiitti 


mwitttmtttMttmttmmttmmttttttmmtttm 


nuiknt  particle  size  is  .haiit  n  radius 

THIS  PARTICLE  COmtIMTES  7.35I264842247Z  SF  TOTRL  ACTIVITY 
FOR  PARTICLES  KTHEEH  .****377  AMR  .*0*231  R  RADIUS 

ACTIVITY/H2  ■  474. 12454276*9  CI/H2 

ACTIVITY/H3  *  .*4496833170329  CI/H3 

AICRAFT  ALTITUDE  «  43M.  METERS 

CABIN  ACTIVITY  *  3.H64219142I9  Cl 

TISSUE  DOSE  TO  EIPOSURE  ♦  9.  HOURS  IS  .4814261512973  RAOS 

DOSE  RATE  *  .27*7638*68933  RA8S-T ISSUE/H8UR 

SKVSHIHE  ROSE  >  .664774*369179  RADS 

SAMI  MFP  «  2*2.216767*766  H  SIMA  I  ■  4224.912631776  H 

mttttmmmmtmmitimmtmtttmimttttttt 


itmtttiiitititiummitiiimitiimtitttittmmtti 


NHIMEVT  PARTICLE  SIZE  IS  .***126  N  RADIUS 

THIS  PARTICLE  C8NTRIDUTES  7.798797*33*62  OF  TOTAL  ACTIVITY 

FOR  PARTICLES  RETUEEM  .****421  AMI  .***231  R  RAIIUS 

ACTIVITY/R2  -  449.22332*7*04  CI/R2 

ACTIVITY/10  •  .*4242744239628  CI/K3 

AICRAFT  ALTITUDE  -  4*H.  METERS 

CABIN  ACTIVITY  -  2.69*I7M7*727  Cl 

TISSUE  I8SE  TO  EIPQOME  *  8.  MORS  IS  .43*7952663343  RADS 

BOSE  RATE  «  .242219476113  RADS-TI99UE/MNR 

SKY9NIRE  ONE  •  .3948*34889961  RADS 

8AMM  HFP  -  191.76*9*93243  i  SIMA  I  -  4224.812631776  H 

tttmittttitmttmitimtttuittitimttmttuimis 


mtttmmtmttittttttattttmmttttttmtmtiimii 


OONINENT  PARTICLE  SIZE  IS  .H**132  H  RADIOS 

THIS  PARTICLE  CONTRIBUTES  R.231M16976942  V  TRTRL  ACTIVITY 

FSR  PARTICLES  IETHEER  .DMD4D4  AM  .1*1249  R  RA9I0S 

ACTIVITY /N2  ■  423.1324333733  CI/H2 

ACTIVITY/ID  >  .*4*13214943601  C1/N3 

AICRAFT  ALTITUK  «  39**.  KTEXS 

MIR  ACTIVITY  •  2.4159433*6329  Cl 

TISSUE  SUE  Tf  E2PM0RE  *  0.  RMS  IS  .306143396*483  IAN 

BASE  RATE  >  .2173697667313  RRM-TISME /HM 

SKYMIHE  f  SE  >  .334! 7386373*8  RAN 

DIMA  RFP  -  181.97913541*2  I  SIMA  I  •  4224. *12631 774  R 

miiiiiittiitiitiitimiimtttmtttittmmtmtmi 
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mtttttittfftftittiftmtmtttttttitttttitiittittttiittt 


•ORIREJIT  PMTT1CUE  SIZE  IS  .99*14  A  RAIIUS 

THIS  PARTICLE  CONTAIWTES  9.7*49499738741  OF  TOT AL  ACTIVITY 

FOR  PARTICLES  9ETKER  ****55  M»  .WMf  H  RMIUS 

ACTIVITY/R2  *  4*2.444474(975  Cl /A? 

ACTIVITY/R3  >  .*3N*2f 24*99233  CI/A3 

AICRAFT  ALTITU9E  -  3 IN.  NETERS 

CMIN  ACTIVITY  »  2.(727342253*9  Cl 

TISSUE  90SE  TO  EXPOSME  ♦  I.  HOMS  IS  .3479*243194*4  RAIS 

ROSE  RATE  *  .1954473459112  AAIS-TISSUE/ ROHR 

SKYSNINE  UK  «  .41*3997143434  RAN 

IRMA  HFF  •  172.1*94125933  9  SUM  I  ■  4224.912431774  A 


ttmtmmttmtmmttmmmtiittttitttmmtmi 


99RIKMT  PARTICLE  SIZE  IS  .9*9141  N  MIMS 

THIS  PARTICLE  CMTRINTES  9.I45231944747X  8F  TtTAL  ACTIVITY 

F9R  PARTICLES  KT9EER  .****591  AM  .9*927  9  AKIM 

ACTIVITY/10  «  392.493577*344  C1/R2 

ACTIVITY/93  •  .93414399991193  CI/RS 

AICMFT  ALT1TNE  -  2599.  STEM 

CAIIR  ACTIVITY  »  1.942133294495  Cl 

TISSK  ME  T9  EIV9ME  ♦  9.  MRS  IS  .3141919995124  9A9S 

ROSE  RATE  •  .1747915251347  RA9S-TISM/IMR 

SKYSNINE  NK  •  .4339341914579  RMS 

SANNA  Wf  -  144.1947427994  R  SIMA  I  •  4224.912431774  R 

ttmtttmttttmttttmtmmimtmtmmmmm 


ttmiiititttimimiimimittimiitiiitttittiiititH 


Mm  PARTICLE  SIZE  IS  .99*154  R  MII9S 

THIS  PARTICLE  C0NTR1IUTES  9.4373*93341*71  9F  TtTAL  ACTIVITY 

FOR  PARTICLES  KTNEEI  .9999453  AM  .99927  R  AMIRS 

ACTIVITY/R2  -  343.9574449993  CI/N2 

ACTIVITY 7R3  -  .13434494939494  CI/R3 

AICMFT  ALTITME  *  299*.  RETERS 

CAIIR  ACTIVITY  •  1. 773439141745  Cl 

TISSK  MK  TO  EIP99M  ♦  I.  MM  IS  .2939995549135  RAM 

I0K  RATE  •  .1997245324493  AA99-T1999E/NMR 

SKYSNIK  IM  •  .3921574799232  RAM 

Mm  RFP  >  154.9979459721  R  SIMR  I  •  4224.912431774  t 

ittmiimtiiitMtMiiiiMAiiiiitsitiismtmiiimiii 
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ttmmitttmtittttitttiittttttttttiitttttmtmttittt 


Ml  ran  PWTICLE  SIZE  IS  .RMI44  H  MI1US 

INIS  PARTICLE  CflHTFIMJTES  lA.f7A47VA44AfZ  OF  TOTAL  ACTIVITY 

for  articles  seneer  .nm7I3  am  .wmi  raiios 

ACTIVITY/R2  •  344.77327M2V5  CI/W2 
ACTIVITY 710  «  .13273141174134  CI/N3 
AICMFT  ALTITME  «  ISM.  KTERS 
CARIM  ACTIVITY  >  I.4MR38V3I42  Cl 

tissue  iose  to  Emmac  ♦  i.  mors  is  .2374123257331  rms 

ROSE  RATE  «  .  1441132414537  RARS-TItHE/MUR 

SKYSNIHE  SHE  >  .  35554272 J74RV  RAIS 

KAMA  MP  «  I4R.4Y41M97M  R  SIAM  I  •  4224.112431774  I 

ttimttttmtttttmtmtmmmimmttiimmmt 


ittmttitmtmitmmtmtmttmiittttmiiiitimi 


IKi£:  "ARTICLE  SIZE  IS  .RM144  N  RARIHS 

THIS  PARTICLE  CSKTRIIUTES  1R.S4M4IV344R1  *  TRIAL  ACTIVITY 

ERR  PARTICLES  SETKEH  .HM 777  AM  .MR2V4  H  RARIHS 

ACTIVITY/IO  >  33R.3SR4V07H3  CI7R2 

ACTIVITY/ID  •  .R3I222I74VC374  CI/K3 

AICRAFT  ALTITME  -  IMS.  REFERS 

CASH  ACTIVITY  -  1.43VI232A4VR2  Cl 

TISSUE  MSE  Tl  EIPRBME  ♦  I.  MRS  IS  .2334374471142  RAM 

ROSE  RATE  '  .  I3I4R2S3Y3R27  RAM-TISSHE/NMR 

SKVSRME  ME  •  .3224UV4M417  RAM 

•AMR  m  •  I4T.3433BAMS  R  SIMA  I  -  4224.  M2431 774  R 

tllltltmilltAtMAIAIAttAAIIIimitllimilllllimitll 


imiitittAAAiMiitttttttimmittttitiMmtitittttmt 


NRIAERT  PARTICLE  SIZE  IS  .IN174  N  RARIHS 

THIS  PARTICLE  C8NTRIMTES  I1.R43S13H3ISR  V  TSTAL  ACTIVITY 

FM  PARTICLES  IETUEEN  .RRRASI4  AM  .RRR323  A  AARIU5 

ACTIVITV/N2  >  3I4.3V2SV1R3I  CI/N2 

ACTIVITY/IO  -  .f2YVflMf 3774V  CI/K3 

AICMFT  ALTITME  »  SIR.  METERS 

CASH  ACTIVITY  -  I.33RB472V4M  Cl 

TISMC  USE  Tl  EtfMME  ♦  I.  RMS  It  .2IJRV7WSA354  RMS 

•ME  MTE  ■  .IIVSSSB44744  AAM-TISSUl/MM 

SKYMRE  ME  -  .2VA2S47333R7  MM 
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Appendix  D 


Microsoft  BASIC  Version 

A  Microsoft  BASIC  version  of  the  computer  code  was  also 
developed.  This  program  performs  the  same  computations  and 
outputs  the  same  values  as  the  FORTRAN  version,  but  was 
developed  for  use  on  a  microcomputer.  The  program  listing 
is  included  in  this  appendix  for  those  who  may  want  to  run 
the  code  on  a  microcomputer. 


79 


II  REN4CAIIN  DOSE  PROWAN  TODCL  21 


21  RERtttOGftM  REVISED  17/11/821 
31  DIN  M(llf),R(lM),Mt(lM> 

41  OPENM/P: ■OUTPUT: DEF  C108U)  *.43431106  (I) 

31  INPin*FFt¥LDlICNfV8ClTRfAGA?*}FFlYlDtACRlVllClTDlAC8 
M  DATA  .48,. 92, I. 29, 1.44, I. 94,2.33,2.71, 3. 17, 3. 45, 3. 83 
71  MTA  4.23,4.44, 5.16,5.44, 5. 94,4.39, 6.17, 7.34, 7.87,8.39 
88  MT8  8.93,9.38,18.1,11.7,11.3,11.9,12.4,13.3,14.1,14.7 
98  DATA  15.4,16.2,17.8,17.9,18.7,19.5,28.5,21.3,22.3,23.5 
188  MTA  24.6,23.7,26.8,21.1,29.3,38.3,31.9,33.2,34.7,34.2 
111  MTA  37.7, 34.3, 41.8, 42.8, 44.6, 44.5, 48.4, 38.5,52.7, 33.1 
121  MTA  57.3, 39.8, 42.5, 63. 3, 48.2, 71.3, 74.3, 77.7, 81. 4,13.2 
138  MTA  89.2, 93.3, 98.8, 183,188, 113, 119, 124, 132, 148 
148  MTA  148, 134, 144,  176,  IN,  281, 213, 231, 249, 27! 

138  MIA  294,323,337,398,438,519,413,762,1833,1873 
148  F08  Il-l  TO  188>8EA8  IlIIHAUlHHlXHlE-fllMIT 
171  T/HCt/VAC!NR-TA/3488iV8L<<244. 337i^4«.2 
188  PtlNTTINE  8F  A8RIVAUNrl*'|M 

190  NCB(44.+6.llt.06(Yl8/l8N)  *9.2831 (IM<Tll/l888)+2. 42)  IA8SU.SB<VIV/I888)*2. 42)  >4384.1 
288  MOOIVU) 

211  Cl«fID17.8899.34tT9.8lt22AIT*2-.885227tV*39.8884J7l9*4) 

228  CKIPU.374-.881197ll*.883434IT*2-.8841t9*3+.888l443ITA4) 

238  PtiNT*a«i  conn  m  iwtid-’inc 
248  8MI-338E«84tnJtFF 
238  B8-2488i8l-i.81323l«8S 
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\  v'V  K.1 V' <■'  <■  ■ 


24#  #MC/3#4.liTCfM2IIC/4#-<2.3l(IC/4#)A2)»TC2l"d.-.5tCtPMICA2)f/(25.*2))iTC,TCttTC2tl.#57J2 
27#  8»<!N(f.7«iH(Tlf/liN)/3)-I.29/44.MLN<VIJ/l«l«)»SJ)A2>mM9.344 

2M  if  «>s.  nw  n-3.  as  th# 

m  si*usr2mi><*.mi/TcnA.3 

3##  FKIPI-.  54  ( (f /SI) A2) )  / <2.3#44tSX) 

31#  NLTHICA/S 

32#  FW  K-l  TO  7 

331  9l*#ilT>3l#i4TVT-#:I»# 

341 

PMITMtPVIffTM,T49(2)  j  *ITEMTI#T|9K(3)  1 1#II0SIW|  SPCI3)  ?  'ST.  (Cl  /H)  '|3PC  (4)  |  'SI.  <N>  'iPKIRTM 

39#  PM  1-1  Tl  IM 

341  HI(I)  42 

37#  Z«C1-C241 

311  THhicZ*#.  IMZ 

3fi  ip  ft-u«ai  m  4i# 

4M  IF  Z>U«9#  m  4*9 

41#  T<2M.I9-#.«raZtP"P«W/2N.19>A5.29#fjllT0  43# 

421  T>214.4SiP>.224IPMnP(-l.3KE-«t<Z-ll#M)) 

431  MID*. ##34i4IP/T}lWMd. IIP /Mdl  )»•#■(!.  43K-#44TAl. 3) /(T*ll#.4)tW*#f/14(  I) 

44#  l^32tllPtT.a#44SI<R(!)A3)/(3MAIimKV*2M:IF  #C>14#  TW#  44# 

43#  KHC/24-2. 334E*#4t  (HCA2)  ♦2.#134€-#4I  <RCA3)  *4.  TIE-###  (1C*4 ) :  MTO  47# 

44#  S-l.2t534*.4t4ia«IS)-.#4U77lia#8(tt)A2»*l.!23S-#3l(aM<S)*3)ttt»l#*«E 

47*  VPHEH(V/»ilC>I.M.I4S-#7/<«mi#UM)i4P-M«K 

4M  Z-Z-V PUT 

If#  TI-TMT 

SC#  IF  TKT#  na  3f# 
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311  IF  Z>RCA»3»SZ  TICI  ill 


32#  IF  Z<AC#*3tSZ  HCN  #2# 

33V  RR(I)*.#!#RAIHFIP)-#.5mHCR)/SZ)*2)/)2.3#443tSZ) 

541  IF  96>f  Tin  371 
33V  IF  MOT*#  1)0  37# 

344  IF  MldKMMt'l)  DO  SV-SZ 
37#  PRINT##,  ItR(nfRR(I),Z 
5##  *TOT»«TOT*RR(I) 

ST#  1*1*1 

###  IF  I-l  DO  L>1 
41#  KIT  I 

#21  PRINTMiPRIRTV#,”!  TIT.  CT.*;9PC(3);*PRRTICLE  •»!«(#>  SPRINT#* 

#31  M-| 
tVFIHTII 
#»  K»*(M(;)/«Tvmi«v 
###  PRINT##, PER, R(J) 

#71  KIT  J 

#a#  iTMiir2)iu.*ii.m)/io*(iitiM»*2r.3tvM. 

Ml  FMVKSt((TI/ST>*2»/t2.3l##tST> 

7W  #2*#TIT/l2.3###t«T) 

71#  R3B#17(2.K6#tSIl 

72#  PRINTMiPRtlTI#, •RCTIfITT /10  ■"|A3)aCMrv”KTIVfrr/l0  -•|A3)aCl/N3* 

73#  PRINT##,  HTT  ALT  ■*|RCR)*N*,'RRiIRRL  TIK  ■•|M;*NNRS' 

74#  PRINT##, (RCFT  KLOCITT  ••|KC|*N/«C' 

73#  IF  RDK*II«lf  TKX  T*2RR  .#l#3tRC#  ELSE  T*21#.#3 

74#  IF  #Cft<all#N  IKR  P-P#m/2N)*3.23  as  P«.224IP#IEIPH.S#2E-#4»<RQHI#f#» 
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77#  K*3.4MC*43lf/T 


7N  CACT*A2tNHKtV#Ct6#> 

m  MAT»CACTI3.7E*l#l.##3l3.l7tl.4E-#4t3.6/(V0ltl##) 
at#  T»«5im#Tniir-.2-»i«»Tf)A-.2> 

II#  PtMTM,*tISTMCE  FMR  MIST  •*jAC*/l###j*tP 
12#  nt!ITM,*C#IIN  ACTIVITY  «-;CACT*“CI“ 

13#  PIIRTM,  *IOSE  MTE  >';#MT;MM#S  TISSUE  PH  MM* 

W  PtirM, 'TISSUE  MSE  T9  EIP06ME  ♦*jTI|*MMS  IS'|Tt|'MK* 

H#  NM-.##3>MT-.#M3MK 

M  IF  l/wr>.|t$x  Tin  MtINTM,’VYMIK  MS  IMCOMTE  Ml  Tl  LAME  MM  NEAR  FKE  NTT 
171  ll>MT9m.?E»l#MM/Mm<FY/«K)<i.tf-ll 
m  SWItMT-1.2) 

IK  PtMTM,*«VMME  MS  -*|9»)*RAMa 
TM  FtlWra^IMM  MFF  •  *|l/MT|a  Na,*MMA  I  >  *|»|*  P 
Til  ADHO-MLT 
V2#  IF  M  Tin  tCMN 
f]#  IF  X-A  Tin  A02S# 

»«•  mn  i 

f»  EM 
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received  between  500  and  12,000  meters  altitude.  The  doses  are 
computed  by  considering  the  cloud  size,  the  aircraft’s  transit 
time,  the  ingestion  rate  of  radioactive  particles,  the  aircraft’s 
distance  to  the  burst  and  the  aircaft’s  altitude.  An  extension 
to  the  computer  program  also  computes  the  dose  received  from  a 
multiple  burst  scenario.  The  results  show,  for  a  single  burst, 
the  total  ionizing  dose  each  aircrew  member  receives  is  approxi¬ 
mately  5  rem.  This  is  for  an  altitude  of  9500  meters  one  hour 
after  cloud  stabilization.  The  multiple  burst  cane,  under  the 
the  same  conditions,  is  approximately  ?04  rem.  both  the  single 
and  multiple  burst  cases  use  a  mission  completion  time  of  8 
hours  after  entering  the  cloud. 
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